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SUMMARY 


This  report  constitutes  a reference  manual  for  the  Underwater  Shock  Analysis  (USA) 
Code,  a computer  program  for  calculation  of  the  transient  response  of  a submerged  struc- 
ture to  a spherical  shock  wave  of  arbitrary  pressure  profile  and  source  location.  The 
code  considers  the  structure  to  be  linear-elastic  and  treats  the  surrounding  fluid  as  an 
Infinite  acoustic  medium.  A discrete-element  (finite-element,  finite-difference)  compu- 
tational model  is  used  for  the  structure,  while  the  computational  model  for  the  fluid  is 
based  upon  the  Doubly  Asymptotic  Approximation. 
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SECTION  1 


INTRODUCTION 

This  report  documents  a computer  program,  the  Underwater  Shock  Analysis  (USA)  Code, 
that  calculates  the  transient  response  of  a submerged  structure  to  a spherical  shock  wave  of 
arbitrary  pressure-profile  and  source  location.  The  structure  is  considered  to  be  linear- 
elastic  and  the  surrounding  fluid  is  treated  as  an  infinite  acoustic  medium.  The  computa- 
tional model  for  the  structure  is  constructed  through  the  use  of  an  auxiliary  discrete-element 
(finite-element,  finite-difference)  code  of  choice  [1,  2],  while  that  for  the  fluid  is  con- 
structed through  the  use  of  the  Doubly  . symptotic  Approximation  (DAA)  [3,  4]. 

1.1  DOUBLY  ASYMPTOTIC  APPROXIMATION 

The  principal  advantage  of  the  DAA  is  that  it  models  the  infinite  acoustic  medium  sur- 
rounding the  structure  as  a membrane  covering  the  wet  surface  of  the  structure.  Hence  fluid 
motion  is  described  merely  in  terms  of  wet-surface  response  variables,  which  are  then  linked 
by  compatibility  relations  to  the  structural  response  variables.  Furthermore,  this  descrip- 
tion is  a simple  matrix  ordinary  differential  equation  with  desirable  computational 
properties . 

The  principal  disadvantage  of  the  DAA  is  that  it  constitutes  an  approximation  to  the 
"exact"  boundary-element  representation  of  the  surrounding  medium  [5,  6J.  The  DAA  does  ap- 
proach exactness  for  both  high-frequency  (early-time)  and  low-frequency  (late-time)  structural 
motions,  however,  and  effects  a smooth  transition  between  the  two  asymptotes.  In  addition, 
it  has  exhibited  satisfactory  accuracy  in  a variety  of  check  calculations  [4,  5,  7].  Hence, 
in  view  of  its  desirable  computational  properties,  the  DAA  is  considered  suitable  for 
engineering  analysis. 

1.2  STAGGERED  SOLUTION  PROCEDURE 

The  governing  matrix  equation  for  structural  response  Is  a second-order  ordinary  dif- 
ferential equation  in  time,  while  that  for  fluid  response  is  a first-order  ordinary  differ- 
ential equation.  Simultaneous  solution  of  these  equations  by  direct  step-by-step  numerical 
Integration,  however,  is  unacceptably  expensive.  Hence  the  USA  Code  utilizes  a staggered 
solution  procedure  [8]  for  step-by-step  solution  of  the  equations  in  time. 

Now  a staggered  solution  procedure  involves  a response  extrapolation  at  each  time  step, 
which  usually  leads  to  numerical  instability  for  time  Increments  exceeding  a critical  value. 
Because  this  critical  value  may  be  unacceptably  small  for  many  computations,  the  governing 
equations  for  fluid  response  have  been  modified  in  such  a way  that  unconditional  stability 
is  achieved.  Thus,  through  avoidance  of  both  direct  simultaneous  solution  and  condition- 
al stablty  constraints,  highly  efficient  computation  is  possible  for  the  greatest  variety 
of  cases. 
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As  an  illustration  of  the  capabilities  of  USA,  a transient  response  calculation  has 
been  performed  for  a 2490  degree-of-f reedom  (DOF)  structural  model  with  a stiffness-matrix 
average  half-bandwidth  of  85  DOF.  The  central-processlng-unlt  (CPU)  time  on  a Univac  1108 
required  for  the  280  time-step  calculation  (with  a single  change  in  time  Increment  during 
the  calculation)  was  28  minutes.  The  corresponding  time  on  a CDC  6600,  on  which  the  code 
also  operates,  would  be  about  10  minutes. 

1.3  INPUT/OUTPUT 

The  USA  Code  requires  three  types  of  input  data  in  order  to  perform  its  function.  First, 
structural  mesh-geometry,  mass-matrix  and  stiffness-matrix  data  must  be  proviaed  by  the 
structural  analysis  code  used  by  the  analyst.  Second,  fluid  mesh-geometry  and  boundary- 
element  data  must  be  furnished.  Finally,  charge  standoff  and  incident  pressure-profile,  as 
well  as  time  integration  specifications  must  be  provided. 

The  code,  in  its  turn,  outputs  structural  displacement  and  velocity  histories  and  fluid 
pressure  histories  for  the  wet  surface.  Response  data  post-processors  furnish  nseudo-veloclty 
shock  spectra,  and  response-history  and  shock-spectrum  plots.  In  addition,  post-processors 
embedded  in  the  structural  analysis  code  may  be  used  to  obtain,  for  example,  stress  and  strain 
response  histories,  as  well  as  response-history  and  stress-contour  plots.  As  currently  con- 
figured, the  USA  Code  can  routinely  handle  problems  with  up  to  2500  structural  and  145  fluid 
DOF  within  a core  allocation  of  65000  decimal  words.  If  primary  core  size  is  increased  to 
100000  decimal  words,  approximately  200  fluid  DOF  may  be  processed  with  a corresponding 
increase  in  the  structural  DOF. 

1.4  SPECIAL  FEATURES 

A number  of  special  features  are  incorporated  in  the  code.  First,  a capability  has  been 
provided  to  handle  a fluid  mesh  on  the  wet  surface  that  is  not  coincident  with  the  surface 
mesh  for  the  structural  model.  This  permits,  for  example,  the  use  of  a refined  structural 
mesh  in  a region  of  high  stress  gradients,  even  though  a relatively  coarse  mesh  is  retained 
for  the  fluid. 

Second,  options  for  variable-step  time  Integration  and  computation  restart  are  furnished. 
The  former  allows  the  use  of  small  time  Increments  during  periods  where  the  response  is  ex- 
pected to  be  varying  rapidly  in  time,  and  the  use  of  large  time  Increments  for  periods 
characterized  by  a slowly  varying  response.  The  latter  permits  the  division  of  a response  com- 
putation Into  segments,  so  that  the  analyst  may  examine  the  results  at  selected  points  along 
the  way.  Such  examination  is  facilitated  by  the  use  of  the  "printer-plot"  routine  that 
augments  the  usual  printout  data  with  response  plots  "drawn"  by  the  printer. 
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Finally,  the  code  Incorporates  fluid  wet-surface  elements  for  both  general  and  beam-iike 
motions  of  the  structure  (see  Appendix  A).  This  feature  is  especially  useful  for  compartment- 
by-rompar tment  analysis  of  <a  submarine.  Such  an  analysis  utilizes  a general-structure 
discrete-element  model  of  a particular  compartment  of  interest,  with  the  remainder  of  the 
submarine  modeled  as  a beam.  Hence  a detailed  analysis  of  an  entire  submarine  may  be  per- 
formed with  several  discrete-element  models  of  moderate  size,  avoiding  the  use  of  a single 
gigantic  model. 

1.5  SPECIAL  NON-FEATURES 

Several  features  of  modest  complexity  have  yet  to  be  incorporated  into  the  USA  Code. 
First,  an  option  for  automatic  time-step  integration  would  free  the  analyst  from  having  to 
select  integration  time  increments  in  accordance  with  his  expectations  regarding  response 
behavior.  Second,  fluid  wet-surface  elements  for  bar-like  motions  of  a structure  would  allow 
proper  treatment  of  the  longitudinal  vibrations  of  a submarine  analyzed  on  a compartment- 
by-compartment  basis.  Finally,  a capability  to  handle  very  large  problems  would  be  useful 
in  those  cases  where  structural  segmentation  is  not  possible. 

Two  important  features  of  greater  complexity  have  yet  to  be  incorporated  into  the  code. 
The  first  is  a treatment  of  hull  cavitation,  which  may  substantially  affect  structural 
response  for  incident  shock  waves  of  short  duration.  The  second  is  the  Inclusion  of  free- 
surface  (bulk  cavitation)  effects,  especially  as  they  influence  the  form  of  the  incident-wave 
excitation.  The  introduction  of  these  two  features  requires  the  treatment  of  highly  non- 
linear phenomena,  which  presents  a stiff  challenge  for  future  work. 
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THEORY 

This  section  tles(  rihes  C(ie  theoretical  foundation  of  tlie  USA  Code.  it  is  constructed 
as  an  overview,  with  coverage  of  details  left  to  referenced  papers  and  reports,  and  to 
appendices . 

2.1  STRUCTUfUU.  RESPONSE  EQUATION 

The  matrix  ordinary  differential  equation  for  the  dynamic  response  of  a linear-elastic 
SI ruci ure  is  [ ] | 


Mx  + Cx  + Kx=f  f2.1) 

'-s  — ~s  - ~s  - — 

where  x is  the  structural  displacement  vector,  .M  , C and  K are  the  structural  mass. 

~s  ~s  ~s 

damping  and  stiffness  matrices,  respectively,  £ is  the  external  force  vector,  and  a dot 

denotes  a temporal  derivative.  Generally,  , C^  and  are  highly  banded,  symmetric 

matrices  of  large  order;  at  present,  the  USA  Code  considers  M to  be  diagonal  and  C to 

■~“S 

be  zero. 

For  excitation  of  a submerged  structure  by  an  acoustic  wave,  is  given  by 

i = -£  ^2.2) 

where  and  p^  are  nodal  pressure  vectors  for  the  wet-surface  fluid  mesh  pertaining  to 

the  (knownl  incident  wave  and  the  (unknown)  scattered  wave,  respectively,  A^.  is  tlie  diag- 
onal area  nuatrlx  associated  with  elements  in  the  fluid  mesh,  and  G is  the  transformation 
matrix  that  relates  the  structural  and  fluid  nodal  surface  forces.  More  will  be  said  about 
G in  the  next  subsection, 

2.2  DAA  EQUATION 

The  Doubly  Asymptotic  Approximation  may  be  written  [1,4) 

Mf  P.S  + f'^AfPs  = (2.3) 

where  u^  is  the  vector  of  scattered-wave  fluid-particle  velocities  normal  to  the 
structure's  wet  surface,  p and  c are  the  density  and  sound  velocity  of  tlie  fluid, 
respectively,  and  Mj.  is  the  symmetric  fluid  mass  matrix  for  the  wet-surface  fluid  mesh 
(see  Appendix  A).  This  matrix  is  produced  by  a boundary-element  treatment  of  Laplace's 
equation  for  the  irrotatlonal  flow  generated  in  an  infinite,  invlscid,  incompressible  fluid 
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by  motions  of  the  structure's  wet  surface;  it  is  fuily  populated  with  non-zero  matrix  ele- 
ments. When  transformed  into  structural  coordinates,  the  fluid  mass  matrix  yields  the  added 
mass  matrix,  which,  when  combined  with  the  structural  mass  matrix,  yields  the  virtual  mass 
m<atrix  for  motions  of  a structure  submerged  in  an  incompressible  fluid  [9]. 

As  mentioned  in  Section  I,  the  approximate  relation  (2 . 3)  is  called  "doubly  asymptotic" 
because  it  approaches  exactness  in  both  the  high-frequency  (early-time)  and  low-frequency 
(late-time)  limits.  For  high-frequency  motions,  >>  |p^|  , so  that  (3)  approaches  the 

relation  = pcu^,  , which  is  the  correct  limit  for  short  acoustic  wavelengths.  For  low- 
frequency  motions,  )^|  <<  |p^|  , so  that  (2 . 3) approaches  the  incompressible-flow  relation 

Aj.  p^  = M|.  Ug  , which  is  the  correct  limit  for  long  acoustic  wavelengths. 

For  excitation  by  an  incident  acoustic  wave,  u^,  is  related  to  structural  response  by 
the  kinem<itlc  compatibility  relation 

= iLx  l-S  (2.^) 

where  the  superscript  "T"  denotes  matrix  transposition.  Equation  2.4  expresses  the  con- 
straint that  normal  fluid-particle  velocity  match  normal  structural  velocity  on  the  wet 
surface  of  the  structure.  The  fact  that  the  transformation  matrix  relating  those  velo- 
cities  is  G follows  from  the  invariance  of  virtual  work  with  respect  to  either  of  the 
wet  surface  coordinate  systems.  Generally,  £ is  a rectangular  matrix  whose  height 
greatly  exceeds  its  width,  inasmuch  as  the  number  of  structural  DOF  usually  exceeds  con- 
siderably the  number  of  fluid  DOF. 

2.3  I.NTERACTIO.N  EQUATIONS 

The  Introduction  of  (2.2)  into  (2.1)  and  (2.4)  into  (2.3)  yields  the  interaction  equa- 
tions 

«S  K + i + Kg  X = - £ Af  ^ 

j (2.3) 

Mj.  pj,  + p c Aj.  Pj.  = P u (G  X - Uj^) 

These  equations  may  be  solved  simultaneously  at  each  time  step  by  the  transfer  of  “ £ Pt 
T 

and  p c M,.  C x to  the  left  sides  of  their  respective  equations.  Such  a procedure  is 
exceedingly  expensive,  however,  because  of  the  large  connectivity  of  the  coefficient  matrix 
Involved.  As  mentioned  in  Section  I,  efficient  computation  is  possible  through  the  appli- 
cation : staggered  solution  procedure  that  is  unconditionally  stable  with  respect  to 

the  i.ioi  c ol  time  Increment. 

The  simplest  implementation  of  tiie  staggered  solution  procedure  recommended  in  [8]  may 
be  effected  as  follows.  is  taken  to  be  diagonal  and,  to  allow  for  the  possibility  that 

.M  mav  have  zero  entries  for  rotational  DOF,  G is  constructed  such  that  only  the  transla- 
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tlonnl  DOF  for  the  structure  couple  with  the  fluid  DOF  [see  (2.4)];  then  the  first  of  (2.5) 

T.. 

may  be  partitioned  to  obtain  0,  2S»  “hlch  may  then  be  introduced  into  the  second  of  (2.5). 
F’remul  1 1 pi  Icatlon  of  tlie  resulting  equation  by  1 ...  1 


Af  Mr 

pc  ~f 


1 

pc 


Af  Ps 


+ (fif  + fi«)  P.S  = 


then  yields 

-(Ss  P-i  + Ay  Uj) 


-Af  i 


(2.6) 


-i  T -i  —1 

where  Dy  = Ay  My  Ay  and  Uj,  = Ay  C £ ~f  symmetric,  and  where  Mj,  is  a diagonal 

matrix  with  each  nonzero  element  given  as  the  reciprocal  of  the  corresponding  nonzero  element 

M^  each  zero  element  mirroring  the  corresponding  zero  element  of  . The  first 

of  (2.5)  and  (2.6)  are  herein  termed  "tlie  augmented  interaction  equations". 

2.4  SPIIKRICAL  INCiDENT  WAVF, 


Each  element  of  the  vectors 


Pi 


and 


for  a spherical  Incident  wave  are  given  by 


(K  . - S V 

t 

' [pc  ^ pRj  Pii^*'^J 


(2.7) 


where  .S  is  the  "charge  standoff",  l.e.,  the  distance  between  the  origin  of  the  incident 
spherical  wave  and  the  nearest  point  on  the  structure's  wet  surface,  Ry  is  the  distance 
from  t tie  origin  of  the  incident  spherical  wave  to  the  iUi  fluid  node  on  the  wet  surface. 
Yy  is  the  cosine  of  the  angle  between  the  vector  corresponding  to  Ry  and  the  wet-surface 
normal  at  the  Itii  fluid  node,  and  Py(t)  is  the  incident-wave  pressure-profile  defined  at 
Ry  = S . For  a shock  wave,  Py(t)  is  discontinuous  at  t = 0 and  the  '"’^^(t)  contain 
singularities. 

In  order  to  remove  shock-wave  singularities  from  liy  in  (2.6),  a modified  pressure 
vector  is  defined  as 


Pm  “ iPl  Ps  (2.«) 

where  P is  a diagonal  matrix  with  direction-cosine  elements  Yy  • The  Introduction  of 
(2.8)  into  (2.6)  and  the  first  of  (2.5),  followed  by  utilization  of  the  second  of  (2.7)  then 
yields  the  modified,  augmented,  interaction  equations 


M X + C X + K X 

~s  - ~s  - ~s  - 


pc 


Af  E„  + (fif  + fis  ) E„ 


-£Ay  [E„+  (i-D  Py] 

- Ay  (c„  X + Kb  if)  " ilEi 


(2.9) 
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in  which  I is  the  identity  m/itrix,  and 


II  = I) 


- (f'  + 

'^S 


- 


- K ' ) r 


(2.10) 


wliere  jR  Is  tiie  diay,onal  matrix  formed  by  the  distances  R^.  Equations  (2.9)  (witl) 

C = 6)  are  the  equations  solved  by  the  USA  Code  to  determine  the  structural  responses 
X and  X,  and  the  wet-surface  pressures  p = (^  - _r)  Ej 
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SECTION  III 


OKCANIX.ATION 

TIk'  USA  Code  has  hefti  written  in  standard  KOK'I'KAN  iV  fnr  use  on  hotli  llnivac  and  (dKi 
computers.  M.icliine  dependency  has  been  isolated  in  one  utility  program  described  below. 
I’rot’.ram  modularity  has  been  strictly  enforced,  wlLli  communication  between  computation.il 
modules  controlled  by  rae.ins  of  a data  management  system. 

Tlie  liasic  structure  of  tlie  code  Is  sliown  In  Fig.  1-1.  Tlie  structural  preprocessor  is 
a separate  code  selected  liy  the  user  t.  provide  the  computational  model  for  the  structure, 
file  skyline  utility  merely  reformats  M and  as  provided  by  the  structural  preprocessor 

for  processing  liy  tlie  USA  Code  (recall  that  is  t.aken  as  zero).  The  ! In  id  mass  preproc- 
essor forms  Aj.  , Mj.  , .and  C , wliile  the  matrix  augmentation  preprocessor  forms  , 

.and  Aj.  c'm^  ' (see  (2.9)).  Tlie  main  processor  is  tlie  time  integrator,  wliicli  forms  £ and 
H and  then  solves  (2.9)  in  step-by-step  f.ashion  using  the  staggered  solution  procedure. 

The  response  postprocessor  provides  tabular  and  graphic  output  for  the  computed  kinematic 
responses  as  well  as  pseudo-velocity  shock  spectra.  Finally,  the  data  manager  controls 
the  flow  of  data  between  processors.  More  detailed  descriptions  of  the  v<arlous  program 
components  follow,  while  information  required  for  utilization  of  the  code  is  contained 
in  Appendices  B througli  E. 

l.i  TliE  DAiA  MANAGER  DMCASI' 

DMCASP  is  a self-contained  utility  module  th.at  functions  as  a manager  of  auxiliary  stor- 
age and  as  the  focal  point  for  all  block  Input/output  activities  (10).  Constituting  tlie 
lowest  level  of  the  NOSTRA  Data  Management  System  (11),  it  carries  out  the  direct  transfer  oi 
data  blocks  between  core  and  peripheral  storage.  (The  terminology  "direct  transfer"  is  used 
here  to  denote  unformatted  and  unbuffered  data  transmission.)  The  basic  auxiliary  storage 
management  operations  embodied  In  DMCASI’  are 

• Activate  storage  device 

• I’ositlon  device 

• Read  data  block  from  device 

• Write  data  block  on  device 

• Deactivate  device 

in  the  USA  Code,  DMCASP  is  operated  as  a stand-alone  1/0  package  that  receives  direc- 
tives directly  from  the  master  processors.  Assembly  language  versions  of  DMCASP  currently 
exist  for  UN  I VAC  1100  EXEC-8,  CDC  SCOPE  l.A,  and  CDC  NOS  operating  systems;  hence  the  USA 
Ciidi-  may  be  used  only  on  these  systems  at  this  time. 
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3.2  THE  STRUCTURAL  PREPROCESSOR 


This  is  a user-provided  code  that  assembles  the  structural  mass  and  stiffness  matrices 
and  generates  information  that  relates  the  internal  and  external  descriptions  of  the  struc- 
tural DOF.  Input  typically  Includes 

• Mesh  geometry 

Coordinate  systems 

Node  locations 

• Element  definitions 

Type 

Connectivity 

• Material  properties 

Mass  density 

Moduli 

• Constraints 

Symmetry  conditions 

Element  external  constraints  * 

Element  internal  constraints 

Fluid  internal  to  the  submerged  structure  must  be  included  In  the  structural  model.  At 
this  time,  USA  treats  only  diagonal  mass  matrices  associated  with  a lumped  mass  representa- 
tion f)f  the  structure,  and  only  single  precision  matrices  can  be  processed. 

1.3  THE  SKYLINE  UTILITY 

This  preprocessor  converts  the  structural  mass  and  stiffness  matrices  generated  by  the 
structural  preprocessor  into  the  Internal  "skyline"  format  required  by  the  USA  time  integra- 
tion processor  |i;;,ll|.  As  there  are  a variety  of  ways  to  store  large,  sparse,  symmetric 
matrices,  virtually  any  structural  preprocessor  that  is  to  be  used  with  the  USA  Code  will 
require  a utility  package  to  change  the  storage  format.  At  this  time,  conversion  utilities 
have  been  written  for  SPAR  [14]  and  NASTRAN  [15].  User  instructions  for  constructing  the 
skyline  utility  for  other  structural  codes  are  given  in  Appendix  F. 

Figure  3.1  shows  2 paths  to  mass  storage  from  the  skyline  utility.  The  SPAR  converter 
uses  UMGASP  for  both  input  and  output,  whereas  the  NASTRAN  converter  uses  unformatted  buf fered 
FORTRAN  commands  for  input  and  DMGASP  for  output. 

Constraints  are  also  handled  differently  in  these  two  utilities.  NASTRAN  provides 
a reduced  stiffness  matrix  which  already  Incorporates  any  prescribed  constraints.  SPAR 
does  not;  however,  USA  has  the  ability  to  apply  constraints  due  to  symmetry  or  attach- 
ment to  ground  during  the  time  integration.  Structural  DOF  that  must  be  set  to  zero  are 
flagged  by  the  skyline  utility  [13]. 
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3.4  THK  Kl.L'II)  MASS  I'KKPkOCESSOR  FI.UMAS 


This  code  constructs  the  fluid  nuiss  nuitrlx  for  a structure  submerged  In  an  Infinite, 
Invlscld,  Incompressible  fluid  by  the  boundary  element  technique  (9).  In  addition.  It  gen- 
erates fluid  mesh  data  and  a set  of  transformation  coeff lcient.s  relating  the  structural  and 
fluid  DOF.  The  Computation  of  these  coefficients  is  based  upon  the  use  of  centroldal  node 
lor  the  fluid  elements  and  the  assumption  of  a bilinear  variation  of  displacement  over  the 
surface  ot  eacli  structur.al  element.  This  assures  that  the  description  of  the  fluid  pres- 
sure forces  in  the  two  mesh  systems  is  statically  equivalent  without  inducing  moments  at 
the  structural  nodes. 

FI.UMAS  contains  a refined  formulation  for  the  fluid  mass  matrix,  which  includes  the 
primary  effects  of  element  curvature.  In  addition,  it  has  the  capability  to  treat  struc- 
tures containing  both  surface-of-revolut Ion  and  general-geometry  components,  as  described 
in  Appendix  A.  ’The  code  can  also  efficiently  construct  the  fluid  mass  matrix  for  a body 
with  two  planes  of  symmetry  by  using  a mesh  which  covers  1/4  of  the  surface.  Symmetric  or 
antl-symraetric  fluid  motions  can  then  be  Imposed  in  the  quadrants  not  covered  by  the  mesh. 
Two-dimensional  "plane-strain"  behavior  of  long  cylinders  can  also  be  simulated  by  another 
branch  In  the  code.  Finally,  a usefu’  diagnostic  tool  contained  within  the  code  is  the 
capability  to  solve  the  fluid-boundary-mode  problem  JM^-u  “ (9). 

'lypical  input  data  for  this  processor  includes 

• Mesh  geometry 

Fluid  Wet-Surface  Mesh 

Structure  Wet-Surface  Mesh 

• Element  definitions 

General  curved  surface 

Surface  of  revolution 

• Material  property 

Mass  density 

• Constraints 

(Quarter  model 

Lung  cylinder 

A detailed  description  of  the  required  input  data  is  given  in  Appendix  B. 

3.5  THE  AUGMENTED  MATRIX  PREPROCESSOR  AUGMAT 

This  preprocessor  accepts  data  from  the  structural  and  fluid  analyzers  to  construct  the 
specific  matrices  required  for  solution  of  the  augmented  Eq.s.  (2.9).  The  output  of  this 
code  includes  not  only  the  required  matrices  in  skyline  form,  but  also  a distillation  of  the 
output  from  both  the  structural  and  fluid  processors.  This  has  been  done  so  that  only  one 
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permanent  file  need  be  referenced  as  Input  to  the  time  Integrator;  this  results  In  Improved 
data  handling  and  core  usage.  Input  to  this  code  Involves  the  following  Information 

• Mass  matrices 

Fluid 

Structure 

• Structural  UOF  correspondence  table 

External  and  Internal  node  descriptions 

Factorization  order 

UOF  reduction  due  to  constraints 

• Fluid  mesh  geometry 

Global  coordinates  of  fluid  nodes 
Direction  cosines  for  nodal  surface  normals 
Areas  of  fluid  elements 

• Fluid/structure  DOF  transformation  coefficients 

• Fluid  material  properties 

Mass  density 
Speed  of  sound 

Although  this  constitutes  a substantial  amount  of  Information,  almost  all  of  It  Is  retrieved 
from  permanent  data  files.  A detailed  discussion  of  the  required  Input  data  is  contained  in 
Appendix  C. 

3.6  THE  TlMi;  INTEGRATION  PROCESSOR  TIMINT 

This  main  processor  constitutes  an  implementation  of  the  unconditionally  stable  staggered 
solution  technique  developed  In  [8).  The  primary  output  is  a permanent  data  file  that  con- 
tains nodal  values  for  structural  displacement,  structural  velocity  and  wet-surface  pressure 
at  every  time  step.  In  addition,  a parallel  file  Is  created  that  retains  restart  information 
at  time  Intervals  dictated  by  the  user.  The  code  has  a variable  time  step  capability  and  can 
treat  a spherical  Incident  wave  of  arbitrary  pressure  profile  and  charge  location.  Finally, 
selected  response  histories  can  be  listed  and  then  displayed  for  Immediate  examination  us- 
ing a "printer-plot"  graphics  package  embedded  both  in  TIMINT  and  in  POSTPR  (see  Sect.  3.7). 

The  computational  strategy  for  the  staggered  solution  procedure  is  embodied  in  the  fol- 
lowing eight  steps,  assuming  the  solution  Is  known  at  time  t: 

(1)  Estimate  the  unknown  structural  restoring  force  x at  t + A t from  the  extra- 
polation of  current  and  past  values 

(2)  Transform  this  extrapolation  into  fluid  node  values  and  form  the  right-hand  side 
of  the  fluid  equation,  which  also  Involves  the  known  Incident  pressure  at  t -t-  A t 
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(3)  Solve  the  fluid  equation  and  obtain  a preliminary  estimate  of  the  total  pressure  | 

vector  at  t + A t ^ 

(4)  Transform  fluid  pressures  into  structural  nodal  forces 

(5)  Solve  the  structural  equation  for  the  displacement  and  velocity  vectors  at  t + A t 

(6)  Transform  the  computed  structural  restoring  forces  at  t + A t into  fluid  tiode 
values  and  reform  the  right  hand  side  of  the  fluid  equation 

(7)  Re-solve  the  fluid  equation  and  obtain  refined  values  for  the  total  pressures 
at  t + A t 

- 

(8)  Save  system  responses 

Steps  1,  3,  and  5 constitute  the  basic  staggered  solution  technique,  while  Steps  2 and  4 are 

required  because  of  the  difference  between  the  fluid  and  structural  surface  meshes.  Ihe  > 

iteration  on  the  fluid  solution  reflected  in  Steps  6 and  7 has  been  added  to  enhance  accuracy. 

Inasmuch  as  the  computation  time  is  overwhelmed  by  the  structural  solution  requirements,  this 
requires  only  a small  increase  in  total  run  time.  The  use  of  a three-point  extrapolation 
method  in  Step  1 also  improves  accuracy,  as  discussed  in  [8]. 

Implicit  integration  algorithms  have  been  used  for  both  the  fluid  and  structural  equa- 
tions. The  former  is  treated  with  the  3-step  Park  method  [16]  while  the  latter  Is  treated 
with  the  "JO"  implementation  of  the  trapezoidal  rule  [17] . 

Typical  input  to  this  processor  includes 

• Charge  characteristics 

Location 

Pressure  profile 

• Time  step  information 

Start  and  finish  times 
Time  increment  values 

• Restart  data 

• Display  directives 

Displacements 
' -locltles 
iressures 

Pseudo-velocity  shock  spectra 

Detailed  user  information  concerning  TIMINT  is  given  in  Appendix  C. 
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J.7  THK  KKSI’ONSi;  I’OK'IPKOCKSSOK  I’OSTI’K 


This  utility  is  responsible  for  the  listing  and  "printer-plot"  graphic  display  of  sel- 
ected system  responses  and  pseudo-velocity  shock  spectra.  Some  of  the  same  capabilities 
are  also  embedded  In  the  TIMINT  processor  for  Immediate  selective  scanning  of  the  output. 
POSTPR,  however.  Is  used  for  more  detailed  examination  of  the  results  at  a later  time.  As 
a complete  display  of  all  structural  and  fluid  DOF  histories  for  even  a moderate  size  prob- 
lem could  run  Into  thousands  of  pages  of  output  care  must  be  exercised  In  the  selection  of 
data  to  be  displayed.  Usage  of  this  code  Is  discussed  In  Appendix  E. 


SECTION  IV 


EXAMPLE  PROBLEMS 

This  section  presents  results  generated  by  the  USA  Code  for  two  idealized  underwater 
shock  problems.  The  structure  studied  in  the  first  problem  is  a hollow  circular  beam  of 
finite  length,  while  that  involved  in  the  second  problem  is  an  infinite,  circular  cylin- 
drical shell.  In  both  problems,  the  structure  is  excited  by  a transverse,  plane  step- 
wave  of  unit  incident  pressure  and  material  properties  are  used  that  correspond  to  a 
steel  shell  immersed  in  water.  The  input  data  are  normalized  so  that  the  density  and 
speed  of  sound  for  the  fluid  both  equal  unity;  hence,  the  density.  Young's  modulus,  and 
Poisson's  ratio  for  the  structural  material  are  taken  as  7.85,  98.125  and  0.3,  respect- 
ively. The  radius  and  wall  thickness  of  the  beam  and  the  cylinder  are  1 and  0.01,  res- 
pectively, while  the  length  of  the  beam  is  9.  In  order  to  assess  the  accuracy  of  the 
computational  results,  selected  response  histories  are  compared  with  those  obtained  by 
other  methods. 

4.1  CIRCULAR  BEAM 

The  response  variable  of  primary  Interest  in  this  problem  is  the  late-time  asymp- 
totic translational  velocity  of  the  structure.  An  analytical  expression  for  this 
quantity  may  be  obtained  from  (2.5)  by  taking  k = v(t),  where  is  the  vector  of 

direction  cosines  relating  the  translational  motions  of  the  structural  nodes  and  the 
direction  of  propagation  of  the  plane  Incident  wave.  (The  elements  of  ^p^  that  pertain 

to  the  rotational  DOF  are,  of  course,  zero.)  The  introduction  of  this  relation  into  the 

T 

first  of  2.5,  followed  by  premultiplication  of  the  resulting  equation  by  Yp^.  then  yields 

T 

where  m = Y„  M..Yt,  ; this  follows  from  the  fact  that  C y_  = K Yn  =0. 
s — Ps  ~f-'-Ps  ~s-‘-Ps  ~s-4’s  — 

After  the  wave  front  of  the  plane  step-wave  has  enveloped  the  structure,  i.e.,  for 

£■1  “ P c Uj  i 

= P c Uj(tl  - t^) 

Hi  ■ 

where  is  the  fluid  particle  velocity  characterizing  the  step-wave,  J[_  is  the  unity  vec- 
tor, the  asterisk  denotes  the  temporal  integral  of  the  quantity  beneath  it,  t^  is  the  vec- 
tor of  incident-wave  arrival  times  for  the  fluid  surface  elements,  and  ]rp  is  the  vector 
of  direction  cosines  relating  the  normals  of  the  fluid  elements  to  propagation  vector  of 
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the  plane  incident  wave.  In  addition,  |£g|  <<  |£g|  for  late-time  motions  (see  Section 
2.2),  so  that  the  second  of  (2.5)  becomes 


i 


£s 


V - Uj),  t»t^ 


(4.3) 


The  introduction  of  this  relation  into  (4.1)  then  yields 


(m^  + m^)  V = - 


(4.4) 


T T T 

where  the  added  mass  Jtljr£  Ypg • But.  from  (4.4),  £ Yp^  = Xp.  so  that  m^  Is  also 

given  as  m^^  = fpilflp- 

With  £^lpg  “ Ip.  the  first  of  (4.2)  yields  Xpg£  = p c U^y^A^l  = 0,  inasmuc  n as 

the  wet  surface  of  the  structure  is  closed.  Hence,  the  right  side  of  (4.4)  vanishes  for 
t > t^,  which  gives  the  expected  result  v = 0.  This  prompts  the  use  of  Integrated  forms 
of  (4.1)  and  (4.3)  (with  quiescent  Initial  conditions),  which  yields.  Instead  of  (4.4), 

(m^  + m^)  V = - 


The  introduction  of  the  second  and  third  of  (4.2)  into  this  equation  then  provides  the 
desired  expression  for  late-time  asymptotic  translational  velocity 


V 


00 


m , + ra 
d a 

ra  + m 
s a 


(4.6) 


T 

where  the  structure's  displaced  mass  m^  may  be  shown  to  be  expressible  as  m^  “ p c Yp 
^j-(tJ_  - Note  that  (4.6)  is  a general  result,  applicable  to  any  wet-surface  geometry. 


Two  different  uniform  mesh  geometries  were  u.sed  to  study  the  circular  beam.  Ten- 
and  twenty-node  models  were  constructed  with  beam  elements  provided  by  the  structural 
analyzer  SPAR  [18].  The  corresponding  fluid  models  contained  9 and  19  elements  of  equal 

size,  with  12  and  24  circumferential  integration  points  (see  Appendix  A).  For  the  beam 

T 

considered,  m^  « 4.439  and  m^  - 28.274;  with  determined  as  m^  ■ X ®esh  geometry 

has  a small  effect  on  the  value  calculated  for  V . It  was  found  that  m « 23.824  for 

oo  a 

the  coarse  mesh  and  m^  • 24.332  for  the  fine  mesh,  which  yield  » 1.843  and  ” 1.828, 
respectively. 

In  the  response  calculations,  a constant  time  step  of  0.1  (20  steps  per  envelopment 
period)  was  used  for  both  models;  the  results  are  shown  in  Figures  4-1  and  4-2.  Velocities 
at  the  ends  of  the  beam  are  higher  than  those  at  the  center  because  the  three-dimensional 
flow  field  at  the  ends  offers  less  resistance  to  the  plane  wave  excitation  than  the  two- 
dimensional  flow  field  at  the  center.  It  is  noted  that  the  responses  of  both  models  are 
similar  although  those  for  the  finer  mesh  appear  to  tend  to  the  rigid  body  asymptotic  velo- 
city more  precisely. 
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4.2  INFINITE  CYI.INDRICAI.  SHELL 


For  this  problem,  a 72-node,  36-eIement  SPAR  model  with  a uniform  circumferential 
mesh  was  constructed.  The  length  of  the  cylindrical  shell  equalled  the  circumferential 
dimension  of  the  square  plate  elements  used  for  the  model;  hence  the  shell  was  one  ele- 
ment long.  Kinematic  constraints  of  zero  axial  displacement  and  no  end  rotation  were 
enforced  through  the  skyline  utility,  as  described  in  Section  3.3.  The  fluid  model  con- 
sists of  36  equally-spaced  elements  around  the  circumference;  the  two-dimensional  nature 
of  the  infinite  shell  geometry  was  simulated  by  exercising  an  option  in  the  fluid  pre- 
processor FLIJMAS  that  adds  flctitiour.  elements  in  the  axial  direction. 

Two  dimensional  n»0,  1,  and  2 modal  response  results  were  generated  by  USA  using  ap- 
propriate Fourier  components  of  the  incident-wave  force  vectors  in  (2.9).  For  comparison, 
liAA  analytical  solutions  were  generated  by  the  method  described  in  [7]  and  fl9).  The 
prlniiiry  response  variables  of  Interest  were  radial  displacement  for  n=0,  radial  and  tan- 
gential velocity  for  n=l,  and  radial  and  tangential  displacement  for  n=2.  A time  step 
of  0.025  was  used  up  to  t”l;  for  t between  1 and  2 this  was  Increased  0.05,  and  for  t 
greater  than  2 a time  step  of  0.1  was  used. 

The  USA  and  corresponding  analytical  results  are  shown  harmonic  by  harmonic  in  Fig- 
ures 4-3  through  4-7.  In  all  cases  the  maximum  errors  made  by  USA  fall  into  the  range  of 
I to  2%. 
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Figure  4-4.  n“l  Radial  Velocity  of  Infinite  Cylinder 
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Figure  4-7.  n-2  Tangential  Displacement  of  Infinite  Cylinder 
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APPEimiX  A 


FLUID  MASS  COMPUTATION  FOR  A SUBMFUGFD 
STKUCTURF  WITH  BOTH  ROTATIONALLV  SYMMETRIC  AND 
GENERAL-GEOMETRY  COMPONENTS 

Computer  calculation  of  the  fluid  mass  matrix  for  a submerged  structure  which  is  composed  of 
both  rotationaily  symmetric  surfaces  and  general  nonsymmetric  surfaces  may  be  accomplished  using 
the  fuily  three-dimensional  boundary  integral  approach  of  Ref.  ( 9 ] . In  such  a case,  however,  the 
use  of  general-surface  elements  for  each  rotationaily  symmetric  surface  is  not  efficient  when  it  is 
known  a priori  that  only  certain  circumferential  Iwrmonics  of  the  fluid  motion  are  of  importance. 

This  may  occur,  for  example,  when  part  of  the  structure  consists  of  a beam  of  circular  cross  section 
for  which  only  the  n = 1 harmonic  is  required.  Then,  the  order  of  the  fluid  mass  matrix  may  be  re- 
duced from  that  for  the  fully  three-dimensional  structure  by  a redefinition  of  the  fluid  degrees -of- 
freerfom . 

For  a general-surface  element  with  a single  centroidal  control  point,  the  associated  degree-of- 
freedom  has  a direct  physical  interpretation,  nameiy,  the  fluid  velocity  normiil  to  the  plane  of  the 
element.  For  rotationaily  symmetric  surfaces,  the  elemental  degrees -of-freedom  become  Fourier 
coefficients  lor  the  circumferontial  harmonics  of  the  surface  velocity  at  prescribed  axial  locations. 
This  gives  rise  to  the  concept  of  the  "surface  of  revolution"  (SOR)  element  as  opposed  to  the  "general- 
surface"  (GEN)  element.  In  this  appendix,  we  describe  a technique  for  the  inclusion  of  both  SOR  and 
GEN  elements  in  the  same  fluid  mass  matrix. 

A direct  way  to  couple  both  SOR  and  GEN  elements  in  a single  matrix  is  first  to  construct  the 
fluid  mass  matrix  for  the  fully  three-dimensional  structure.  Then,  a transformation  can  be  applied 
which  Is  associated  with  a hybrid  surface  velocity  vector  representing  both  physical  and  generalized 
c(X)rdlnates.  To  Implement  this  two-stage  procedure,  we  cover  the  entire  surface  of  the  structure 
with  GEN  elements  and,  in  particular,  require  equal  spacing  of  such  elements  around  the  circumfer- 
ence of  the  surfaces -of -revolution.  The  number  of  such  elements  in  this  latter  category  is  taken  to 
be  the  same  at  all  jixlal  stations  and  is  denoted  by  N . A set  of  N circumferential  GEN  elements 
at  any  axial  station  is  then  regarded  as  a single  SOR  element  having  as  many  degrees-of-freejdom 
as  the  number  of  Fourier  coefficients  which  are  chosen  to  describe  the  global  fluid  motion.  The  set 
of  control  points  associated  with  the  circumferential  GEN  elements  will  be  referred  to  as  integration 
points  for  their  associated  SOR  elements. 

We  now  assume  that  the  fluid  velocity  vector  for  the  full  three-dimensioniil  surface  discreti- 
zation V can  be  approximated  in  the  mean  by  the  hybrid  fluid  velocity  vector  y*  by  the 
transformation 
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where  the  superscript  G refers  to  the  GEN  elements  located  on  the  general  geometry  portions  of  the 
structure  and  the  subscript  S refers  to  the  sets  of  GEN  elements  which  make  up  the  individual  SOR 
elements;  the  vector  v contains  the  Fourier  velocity  coefficients  for  these  SOR  elements.  For  con- 
venience, we  have  partitioned  the  above  arrays  so  that  the  entries  pertaining  to  GEN  elements  appear 
first  and  will  continue  this  convention  in  the  following.  The  matrix  F is  rectangular  and  has  the  form 
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Here,  M is  the  maximum  number  of  circumferential  harmonics  and  L is  the  total  number  of  SOR 

elements  while  the  vectors  and  for  the  ! SOR  element  and  the  m circumferential 
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where  i is  the  circumferential  integration  point  index.  If  the  transformation  (A.  1)  is  now  applied 
to  the  fluid  velocity  vector,  then  invariance  of  the  total  kinetic  energy  of  the  fluid  provides  the  direct 
means  for  reduction  of  the  mass  matrix  as 

M*  = M D (A. 6) 

whore  the  transcipt  T denotes  matrix  transposition  and  M is  the  previously  computed  fluid  mass 
matrix  for  the  fully  three-dimensional  structure. 


Application  of  the  preceding  procedure  can  involve  a large  expenditure  of  computational  effort, 
particularly  if  the  number  of  .SOR  elements  and  circumferential  integration  points  is  large.  In  such  a 
case,  matrix  partitioning  and  out-of-core  operations  may  be  reejuired.  Even  when  large  matrices  can 
be  handled  by  such  core  mamigement  techniques,  the  computations  include  the  factorization  of  a large, 
full,  unsymmetric  matrix  and  subsequent  solution  operations  just  to  form  the  mass  matrix  before  re- 
fluction.  These  considerations  have  prompted  the  development  of  an  alternate  but  approximate  ap- 
proach that  uses  intermediately  formed  matrices  of  the  same  order  as  the  final  reduced  mass  matrix 
with  a corresponding  saving  of  computation  time. 


Calculation  of  the  fluid,  mass  matrix  based  upon  the  simple  source  formulation  involves  three 
state  variables,  the  source  strength  a • the  velocity  potential  0 and  the  normal  derivative  of  the 
velocity  potential  , which  are  related  by  the  matrix  equations  f 9 ] 
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Here,  B anrl  C are  full,  unsymmetric,  square  matrices  while  , and  a are  column  vectors. 

We  take  the  order  of  this  system  to  be  the  number  of  discrete  fluid  degrees -of-freedom  for  the  full 
three-dimensional  surface  discretization.  The  actual  coefficients  of  B and  C and  the  precise 
manner  in  which  they  are  combined  to  form  the  fluid  mass  matrix  are  unimportant  to  the  development 
here  and  the  interested  reader  is  referred  to  | 9 ] for  this  exposition.  However,  we  will  use  cer- 
tain aspects  of  the  connectivity  of  B and  C later. 

A fundamental  assumption  made  in  ( 9 ] is  that  the  source  strength  a is  constant  over  each 
element.  Here,  we  make  the  same  assumption  for  all  elements,  anfl  further  stipulate  that  the  vari- 
ation of  a awl  0 as  well  as  (the  negative  of  the  fluid  velocity  vector)  around  the  circumference 
of  each  SOR  element  can  be  expressed  by  the  same  finite  Fourier  series  embedded  in  the  transforma- 
tion in  (A.  1).  If  there  are  only  SOR  elements  In  the  structure,  it  can  be  shown  that  this  relationship 
is  exact,  harmonic  by  harmonic,  however,  the  presence  of  general -geometry  structure  components 
in  the  flow  field  perturbs  this  exactness.  We  then  have 
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Wc  an-  now  able  to  construct  a sot  ol  f^juations  in  the  roducod  systom  similar  to  Kq.  (A. 7)  by 
multipiyinK  Kq.  (A.'J)  by  the  lo/t  inverse  of  O , which  is  defim^d  as  ( I '< ) . 
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It  then  follows  that  the  counterparts  ol  H ami  C in  the  reducefi  system  are 
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is  easily  obtaimsl  from  the  ortho){onality  and  normalization  conditions  for  the  linitf;  Kourier 
series,  i.e. , I 20  ] 
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IIsinK  K'lS . (A. 2),  (A. 4),  (A.ri),  and  (A.  12),  wo  find  that 
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whirre  G is  a rlia);onal  matrix  Kiven  by 
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w<;  IKIW  show  that  lh<;  formation  of  B anfl  C and  tho  transformation  inilicatnd  in  Bq.  (A.  1) 
l an  rarrifrd  out  si mullanoousiy  to  ({tvo  ho  matrioos  B*  ami  C*  flirootly.  Consistent  with  the 
[)a rtitionin;'  of  Kqs.  (A. 2)  and  tA..'i),  B may  ho  expressed  as 
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where  the  first  sup(!rscript  rcders  to  the  element  type  (flKN  or  Sf)B)  for  which  the  potential  is  beinK 
<!valuated  ami  the  second  superscript  indicates  the  element  tyjK;  whose  source  strentith  contributes 
to  that  potential  ( gee  the  first  of  Kf) . (A. 7)].  Now,  usin^Kfls.  (A..'i),  (A. 10),  (A. 11),  (A.Bi),  and 
(A . 1 5) , we  obtain 
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A similar  result  can  b»r  written  for  C*  ; howfwer,  as  the  development  to  follow  is  the  same  for  both 
B*  ami  C*,  only  that  for  B*  is  presenbtd. 

liquation  (A.  10)  shows  that  th«-  relationship  bjitween  CIKN  fdements  is  unalt(!red  fjy  the  trans- 
formation, but  that  the  matric<'s  couplini?  fi FN  and  .SOB  elements  require  matrix  multiplication.  For 
the  matrix  F , w<:  have 
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wher»r  the  subscript  I identifies  a partirrular  CKN  element,  the  subscript  j identifies  a particular 

SOIt  element,  ami  the  imiex  J'  irientifies  one  of  the  N lnte>{ratlon  fioints  around  the  circumference 
th  i' 

of  the  J SOK  element.  Hence,  the  bjj  can  be  computed  se<|uentially  over  the  Index  j*  and  the 

summation  [Mtrformed  conjunctively.  The  same  technique,  of  course,  can  be  applied  to  the  matrix 
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Construction  of  the  matrix  relating  SOR  elements  is  complicated  by  the  fact  that  double  sum- 
mations are  involved,  f’rom  Eq.  (A.  16),  we  have 
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where  the  subscripts  i and  j identify  SOR  elements  and  the  indices  i'  and  j'  identify  integration 
points  on  the  i^*'  and  / SOR  elements,  respectively.  It  will  now  be  shown  that  (A. 20)  can 
be  simplified  to  an  equation  involving  a single  summation  by  examination  of  the  connectivity  of  B‘  ‘ . 

As  discussed  in  | 9 ] , the  matrix  elements  of  B and  C are,  except  for  a multiplicative  con- 
stant, functions  only  of  the  geometry  of  the  finite  element  mesh.  Now,  it  is  clear  from  Figure  A-1 
that  the  relative  geometry  of  |;)oints  A/D  is  identical  to  that  of  points  B/E,  that  the  relative  geome- 
try of  points  A/E  is  identical  to  that  of  B/F,  etc.  Hence,  the  N x N submatrix  that  relates  the 
source  strengths  at  the  integration  points  on  segment  j to  the  values  of  the  potential  at  the  inte- 
gration points  on  segment  i must  have  identical  terms  along  every  diagomil.  Moreoever,  the  vari- 
ation in  any  row  (or  column)  is  that  of  an  even  function  because  of  the  rotationally  symmetric  geom- 
etry, i.e. , the  relationship  of  C/D  is  identical  to  that  of  C/F,  etc.  On  the  basis  of  these  facts, 

Eq.  (A. 20)  may  be  considerably  simplified. 

Consider,  for  example,  the  expansion 
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with  the  condition  that 
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anrl  the  understanding  that  if  a subscript  on  the  left  becomes  greater  than  N it  is  replaced  by 
k 1 - N,  etc.  Now,  the  double  summation  in  Eq.  (A. 21)  is  ordered  by  rows  and  columns;  however, 
we  can  rearrange  it  Into  N terms,  each  of  which  represents  summation  along  either  the  main 
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diagonal  or  an  extended  diagonal  . Then,  taking  advantage  of  Eq.  (A. 22),  we  obtain 
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Next,  we  reorganize  the  general  term  that  appears  in  Eq.  (A. 23)  as 
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II  m ?!  n , it  follows  from  Eq.  (A.  12)  that  the  general  term  is  zero  and,  hence,  the  sum  S 
equal  to  zero.  If  m = n , then  Eqs.  (A. 12)  and  (A. 24)  reduce  Eq.  (A. 23)  to 
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In  a similar  manner,  the  other  possible  combinations  of  trigonometric  functions  in 
Eq.  (A. 20)  may  be  simplified;  such  simplification  yields,  for  m ?:  n , 
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i he  latter  consists  of  a fliagonal  in  the  upper  triangle  plus  its  complement  in  the  lower  triangle. 
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while  for  m = n in  the  cos/cos  and  sin/sin  combinations  in  Eq.  (A. 20), 
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and,  for  m = n in  the  sin/cos  and  cos/sia  combinations  in  Eq.  (A. 20), 
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Finally,  as  it  has  been  demonstrated  that  the  must  be  even  functions  around  the  circumference 
of  any  SOR  segment,  the  contribution  of  Eq.  (A. 28)  to  the  final  results  must  vanish,  leaving  only 
Eq.  (A, 27).  This,  of  course,  also  means  that  the  computation  required  to  form  Eq.  (A. 27)  can  es- 
sentially be  halved. 


The  results  of  the  preceding  development  may  be  summarized  as  follows: 

• The  coupling  that  occurs  between  SOR  segments  occurs,  harmonic  by  harmonic,  without 
coupling  of  different  harmonics 

• The  coupling  is  the  same  for  both  the  sine  and  cosine  functions,  and  is  given  by  Eq.  (A. 27) 

• The  sine  and  cosine  functions  of  a particular  harmonic  do  not  couple 

Thus  far,  we  have  substantially  reduce  the  computations  required  to  form  the  B*  and  C*  matrices. 
A few  steps  remain,  however,  in  order  to  reach  the  final  objective,  viz. , efficient  formation  of  the 
fluid  mass  matrix  for  the  coupled  structure. 

As  discussed  in  [ 9 ] , B*  and  C*  are  used  to  form  the  fluid  kinetic  energy  expression, 
which  is  derivable  from  an  appropriate  surface  integral.  In  the  discretized,  fully  three-dimensional 
system,  this  requires  a diagnonal  matrix  dA  whose  elements  are  simply  the  incremental  areas  of 
each  surface  element.  The  kinetic  energy  T is  then 


T = -i  P J dA  0 (A. 29) 

s 

where  P is  the  fluid  density  and  the  integration  extends  over  the  fully  three-dimensional  surface  S . 
We  may  write  this  expression  in  the  reduced  system  by  simply  substituting  from  Eq.  (A. 8),  with  the 
result  that 

cM*  = d^  D (A. 30) 


Carrying  out  the  matrix  multiplications,  we  obtain 


dA* 


O 

T <? 

O F d^  F 


(A. 31) 
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Kor  each  sulenatrix  appropriate  to  a particular  SOR  element  and  harmonic  th.it  appears  in 
T S 

K (1  A 1'',  the  (expansion  nivc-s 


d a*  = 
II 


N 


d a 


n 


(A.:w) 


where  da.,  is  the  incremental  area  for  a sinf^le  intejfration  element  around  the  circumference  of 
the  setiment  in  the  full  three  dimensional  discretization. 


The  precedinj;  results  can  now  be  applied  to  any  of  the  computational  schemes  described  in 
( 9 1 . As  is  noted  there,  however,  one  particular  scheme,  labeled  "12,"  has  been  found  to  bo  sur- 
pri.sinKl,y  accurate  in  relation  to  its  computational  demands  and  hence  has  been  chosen  here  as  the 
"baseline"  methixl.  In  this  case,  the  incremental  areas  appearintj  in  Kqs.  (A. 29)  throu(?h  (A. .'12)  may 
be  replaced  by  the  total  area  of  the  element,  as  the  "12"  scheme  essentially  assumes  that  0 and 
as  well  as  o , are  constant  over  each  surface  element. 

A modest  study  has  been  performed  on  a simple  model  in  which  the  reduced  mass  matrix  has 
been  obtained  directly  from  !■>(.  (A. 6)  and  by  the  approximate  methfxl  outlined  above.  This  model 
consists  of  two  cylinrlrical  sh.dl  segments  with  a radial  "fin"  joined  to  one  segment  on  the  lateral  .sur- 
face a.s  slxiwn  in  Figurr;  A-2;  only  the  n=l  Ix'am  modes  were  kept  in  the  reduced  matrix.  Com- 
parisons have  been  made  of  the  eigenvalues  of  the  "fluid  boundary  mode"  problem  defined  by 


M»  X = A A*  X (A.. 3:5) 

(In  ( 9 I the  \ are  identified  with  moflal  kinetic  energy  comixmehts.) 

Typical  results  from  the  study  are  shown  in  the  following  table: 

A -Direct  A-Approxlmatc  Mode 


0 . 5967 1 

0.56960 

S 

0.478;!8 

0.47826 

A 

0.29624 

0.2.5955 

S 

0.27412 

0.27108 

A 

0.21624 

0.21622 

A 

0.196.34 

0.175.54 

S 

where  S denotes  a symmetric  or  breathing  mode  relative  to  the  two  sides  of  the  "fin"  and  A denotes 
an  antisymmetric  or  rigid-body  "fin"  mode.  We  see  that  the  eigenvalues  for  the  rigid  modes  are  in 
goofi  agreement;  however,  the  eigenvalues  for  the  symmetric  modes  exhibit  discrepancies  of  .I  to 
14  percent.  These  discrepancies  may  be  due  to  the  retention  of  only  the  n = 1 modes  in  the  re- 
fluced  matrix;  only  further  evaluation  under  more  general  conditions  can  clarify  this  point.  At  this 
time,  however,  it  appears  that  the  technique  Introduced  here  Is  a viable  one,  particularly  for  very 
large  problems  whose  size  precludes  application  of  the  direct  method. 
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API'KNDIX  B 

(ISHR  rNFOKMA'CrON  f’OK  TIIIC  f'M/ri)  r'KKI'ROCKSSOR  KI.IJMAS 

rfi(*  frjllowlnj'  lnf'IiKi<*K  nifxll  f i r.'i t lonH  which  were  made*  under  Hcparatc  contract  to 
the  Naval  Surface  Weapons  Center,  Contract  Number  Nb0921 -77-C-On 2 . 

This  appendix  Includes  a copy  of  the  users  manual,  and  a sample  Input  deck  and 
Huhse<|uent  output  for  the  Infinite  cylindrical  shell  problem  presented  In  Section  A. 
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The  following  discussion  is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  Included  here. 

First,  the  terms  appearing  under  "Fluid  Mesh  Geometric  Arrays"  are  defined  as: 

- number  of  corner  points  for  a particular  fluid  element 

- global  cartesian  coordinates  of  the  fluid  element  centroidal 
control  point 

- components  of  the  outward  unit  normal  vector  for  the  fluid  element 

- number  of  struc'  jral  node  points  that  are  coupled  to  a partlcul.ar 
fluid  eleinent  fir  the  purpose  of  force  application 

A00,A20, A1 1 ,A02  - area  and  momentsand  product  of  inertia  of  fluid  element.  Used 
Internally  for  construction  of  the  fluid  mass  matrix  and  of  the 
fluid-structure  transformation  coefficients  for  general  elements. 

BII.CII  - diagonal  terms  of  B and  C matrices  used  for  construction  of  fluid 

mass  matrix  (see  (9|) 

When  SOR  elements  are  included  in  tlie  fluid  mesh  the  following  terms  will  appear  in 
the  output: 

RAD  - radius  of  fluid  element  control  point  from  axis  of  revolution 

NCIR  - number  of  integration  points  used  in  cirdumf erential  direction 

NAXI  - see  NAXIS  in  user  manual 

NFRK  - see  NFREE  in  user  manual 


NCOR 

X.Y.X 

NX,NY,NZ 

NTRA 


Local  Fluid-Structure  Transformation  Coefficients  appear  next.  This  is  a summary 
that  Indicates  which  structural  nodes  couple  with  a particular  fluid  control  point  and 
the  weighting  factor  for  each.  The  weighting  factors  must  always  sum  to  unity  for  any 
fluid  control  point. 

The  eigenvalues  and  eigenvectors  that  follow  the  listing  of  the  added  mass  matrix 
correspond  to  the  "Fluid  Boundary  Mode"  problem  [9].  For  the  infinite  cylindrical  shell 
problem  presented  here,  the  exact  eigenvalues  should  behave  as  1/n  with  corresponding 
modes  cos  nO  and  sin  nO.  The  first  eigenvalue  listed,  .11838+04,  Is  an  approximation 
to  “>  for  n=0  and  it  can  be  seen  that  the  subsequent  eigenvalues  are  relatively  well 
behaved . 


B-13 


Q 

Z 

-J  fO 


Z O I 
lA. 

z 


•^8 

Ul 

O 


< 

^ O • I 
%h  ^ 


> 

u 


<p 

n 


H>  O • »-  I 


z 

< 

ft, 


</)0  0«nM  ID 

i/>  ON  r>  < 
S 


U 


^c»r»«iDiONCDoe^Mn« 


o 

X 


B-14 


i 


PLANE  strain  simulation  OF  INFINITE  CYLINDER 


NJ 

7 


oooooooor  ooo 
o o o rj  t>  o o c>  o r-  o o 

o o o o f j o o o o o o o 

oooooooooooo 
oooooooor.  ooo 
O O O O O o C.>  O O O f.j  o 

O O c:)  O O O O c>  o o o o 

oooooooooooo 


ooooooooooo 
o o o o o o o o o o c> 
ooooooooooo 
ooooooooooo 
CjOOOOC  ooooo 
o c»  o o o o o o o o o 
ooooooooooo 
ooooooooooo 


oooooooooo 
o o o o o o o o o o 
o o o o o o o o 
orjoooooooo 
oooooooooo 
O O o o o o o o o 
oooooooooo 
oooooooooo 


o 

O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

•- 

0 

0 

0 

0 0 

o 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 0 

+ 

4 

f 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 4 

o 

CD 

n 

'll 

V- 

-f 

0 

f'4 

m 

0 

'X> 

n 

r^ 

t'- 

«7 

•7 

t'* 

0 

o> 

'7 

«- 

•D 

t'- 

Oi 

f~ 

•- 

'7 

0 

r- 

'7 

m 

fT»  «T 

o 

w- 

0i 

'D 

»7 

'7 

( * 

0 

1 . 

fD 

•■7 

»7 

'V 

0 

•- 

f4 

»- 

')> 

I/I 

ro 

rn 

'P 

f- 

0 

f-~ 

'0 

n 

•7  10 

o 

rn 

O 

n> 

I' 

'7 

in 

0 

r> 

rx 

in 

'I 

J- 

0 

0 

rti 

in 

(T) 

0 

n> 

1 - 

't 

m 

CM 

t * 

0 

r • 

in 

*7  4- 

o 

'■4 

(;» 

Ct) 

n> 

C) 

{ i 

0 

rj 

'7 

ro 

0 

•t 

'7 

n 

'T 

' 4 

0. 

r> 

< j 

0 

<7> 

rN 

»7  (O 

( ) 

'D 

r) 

Oi 

t - 

0 

T J 

'D 

’J) 

0 

'I) 

0 

T) 

• - 

r> 

0 

•F| 

TD 

n- 

t • 

r> 

' j 

■l> 

'll 

(> 

f7> 

ID 

0 

0 4' 

o 

•0 

iD 

•7 

0 

’7 

r7) 

'n 

'n 

1 - 

n 

<4 

n« 

ri 

•D 

n> 

'7 

0 

'7 

0> 

*0 

'0  01 

o 

t ~ 

-7 

n> 

'T 

•-U 

’X> 

'n 

0 

cn 

n 

'0 

x» 

'7 

0 

'7 

1 - 

1 • 

'7 

n> 

'7 

•0 

<D 

a 

0 

00 

rn 

10 

•O  *7 

o 

*" 

'7 

1- 

'J) 

o» 

0‘ 

o> 

a* 

rn 

f- 

>n 

a> 

r*’/ 

f’) 

'7 

(X; 

( - 

(U 

C' 

<}i 

(Ji 

ni 

(X> 

to 

I I i I I I I I I I I I I I 


0 

0 

0 

0 

0 

0 

0 

0 

4* 

0 

0 

0 

0 

0 

0 

0 

0 

rj 

0 

r> 

0 

0 

0 

t- 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ID 

0 

0 

0 

Cj 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 0 

t 

A 

• 

t 

1 

4 

t 

* 

« 

1 

4 

4 

4 

4 4 

in 

ri 

1/. 

r4 

'D 

(7- 

in 

'D 

ri 

1 - 

Ty 

ri 

cr. 

<■4 

If. 

r > 

■D 

»♦ 

cr. 

fTi 

l/i 

(0 

m 

'7 

r> 

'■ 

<N 

IP  f4 

0 

f - 

fp 

*7 

n 

•D 

0 

•- 

in 

•- 

•- 

n. 

1/. 

♦7 

'n 

ID 

1 - 

0 

iP 

'7 

* J 

X’ 

•- 

r4 

Ifi 

*- 

0 

0> 

iP  n 

0 

1 - 

4 

in 

•7 

1 . 

Tj 

Cj 

m 

r < 

'O 

0 

”7, 

I • 

'7 

If. 

r Y 

1 - 

» p 

' • 

rj 

./ 

1 • 

r^ 

c* 

7 

1 - 

0 

p. 

1'  '7 

0 

n> 

ri 

'7 

X 

( 1 

f , 

■7 

'O 

'■4 

Tl. 

,*■/ 

•T 

T i 

'* 

ry 

y 

0 

pi 

r , 

T 

T 

0 

ri 

•7 

P 

'7 

r4 

P‘ 

'O  *7 

0 

*7) 

■p 

0 

0 

0 

0 

'0 

“ , 

0 

7 

1 - 

0 

13 

0 

F , 

f .. 

Cj 

rj 

t 

p. 

t - 0 

0 

r7» 

'0 

(P 

rsi 

0 

CN 

m 

ro 

m 

IN 

n> 

CM 

'n 

Ip 

n* 

'7 

ry 

r/i 

Ip 

Ip 

IN 

0 

IN 

m 

IP 

IN 

m 

04  (P 

0 

CO 

CD 

'D 

'7 

0 

n 

r - 

m 

{- 

cn 

•7 

c 

'•C 

0 

CO 

rp 

X 

X 

•7 

0 

*7 

t- 

'7 

o> 

'7  IP 

O) 

a> 

00 

7'- 

IP 

in 

m 

“* 

m 

't 

'P 

r- 

CO 

pi 

cn 

0. 

cn 

CO 

t- 

(P 

m 

fO 

•7 

m 

1 

Ip  O* 

I I I I I I I i I I I I I i I i I I 


o o o o o o 
o o o o o o 
fsl  o o o o o o 
o o rj  o o o 
o o o o o o 
O O O O (J  o 
o o o o o o 
o o o o o o 


oooooooooo 

oooooooooo 

oooooooooo 

oooooooooo 

oooooooooo 

cjooooooooo 

oooooooooo 

oooooooooo 


oooooooo 

oooooooo 

oooooooo 

ooooooc-o 

oooooooo 

oooooooo 

oooooooo 

oooooooo 


oooooooo 

oooooooo 

oooooooo 

ooOoooooo 

ooooooooo 

000000f>00 

ooooooooo 

ooooooooo 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0* 

0 

0 

0 

0 

0 

0 

0 

0 

V- 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 0 

4 

4 

4 

4 

4 

4 

4 

4 

4- 

4 

4 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 ♦ 

0 

CO 

»7 

pi 

-p 

0 

ri 

in 

0 

'P 

'n 

r-g 

h* 

n 

't 

r>. 

0 

C7* 

't 

«- 

•D 

r- 

!?• 

1 * 

•- 

-p 

0 

h- 

'P 

m 

cn  •T 

C> 

«- 

I’ 

P 

n 

'P 

• ~ 

0 

I » 

'D 

•7 

-7 

•D 

0 

•- 

1 . 

0 

r- 

»- 

a. 

in 

m 

rn 

Ip 

f- 

0 

I’. 

CD 

•p 

•P  '0 

0 

T' 

0 

P“ 

f . 

't 

m 

IN 

r*- 

0 

: ' 

ri 

ifi 

-7 

0 

0 

T) 

.*) 

T> 

0 

IP 

•p 

m 

04 

I-' 

0 

f'* 

04 

m 

'P 

0 

t 

04 

m 

rt) 

n 

fN 

m 

0 

0 

0 

n\ 

IN 

n 

on 

0 

04 

♦f 

ID 

'I 

04 

D) 

03 

04 

17' 

0 

0 

0 

<n 

CM 

•t  CO 

0 

ID 

0 

0 

I- 

0 

0 

Ip 

CD 

0 

m 

CP 

0 

0 

r-- 

0 

0 

CD 

ID 

Ip 

0 

0> 

0 

0 

Ip 

CD 

0 

CD 

CP 

0 

0 0- 

0 

m 

IN 

cn 

04 

IP 

CP 

<7i 

n 

0 

•p 

01 

Ip 

CP 

04 

0 

r-g 

rn 

t'- 

rn 

04 

cn 

04 

CP 

CP 

<n 

♦P 

0 

'P 

m 

<0 

CP  CM 

0 

I-'. 

n 

0i 

1 

Ip 

IP 

rn 

CO 

0 

CO 

rn 

CP 

P 

•p 

0 

-p 

t-- 

•- 

f" 

-P 

o> 

•P 

p 

CP 

m 

CO 

0 

CO 

m 

'D 

Ip  ^ 

0 

n 

't 

ID 

I'. 

CO 

0» 

cn 

•- 

o> 

0i 

CO 

1-- 

CP 

m 

rn 

rn 

rn 

rf 

CP 

t- 

OJ 

O) 

0) 

o> 

cn 

03 

7-  IP 

I I I I I I I I I I I I I I 


0 

0 

0 

0 

0 

0 

0 

0 

f'. 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0. 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 0 

4 

4 

4- 

4 

4 

4- 

4 

4 

4 

1 

4 

4 

4 

4 

4- 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 4 

0 

in 

04 

«- 

X 

04 

X 

cn 

X 

X 

04 

t'. 

ro 

04 

ro 

IN 

X 

0 

X 

ri 

rp 

ro 

ro 

X 

ro 

X 

'p 

0 

cn 

OM 

X OM 

0 

t- 

‘P 

n 

X 

0 

V- 

If) 

— 

rn 

in 

rn 

•0 

0 

0 

X 

•p 

n 

X 

0 

04 

X 

«- 

0 

cn 

X rt 

0 

r^ 

04 

in 

•p 

f' 

0 

0 

ro 

f4 

ro 

0 

a. 

f' 

If) 

04 

1' 

0 

1' 

f4 

X 

n 

f' 

0 

0 

ro 

t 

ro 

0 

cn 

h*  «P 

0 

0 

0» 

04 

^p 

D 

0 

rg 

•p 

n 

n 

04 

Oi 

P) 

'T 

04 

<3, 

0 

0 

0 

in 

Ov 

n 

ro 

0 

04 

•p 

CJI 

'p 

04 

in 

CD  ^ 

0 

CO 

X 

0 

0 

0 

0 

X 

m 

X 

0 

0) 

f* 

0 

c> 

X 

(T. 

0 

cr 

X 

0 

0 

r - 

0 

0 

X 

1- 

X 

0 

in 

h-  0 

0 

•p 

m 

X 

X 

04 

0 

CM 

rn 

CD 

rn 

04 

<n 

04 

X 

X 

cn 

n 

0 

*7 

cn 

X 

X 

04 

0 

04 

rn 

X 

rn 

04 

cn 

CM  X 

0 

CO 

n 

X 

CO 

•P 

0 

*P 

X 

'P 

O) 

n 

X 

•0 

rn 

ro 

0 

ro 

rn 

X 

p 

•P 

0 

•P 

o* 

0- 

r^* 

«P 

cn 

^ X 

•- 

(7) 

CJ* 

X 

t'* 

X 

X 

rn 

•- 

•- 

rn 

-p 

X 

o* 

CO 

cn 

cn 

<n 

cn 

00 

0- 

X 

X 

m 

•- 

•- 

n 

X 

I I I I t I I I I I I I I I I I I I 


M-r^r>cin(DNCD<y*o  — f'<^*Tm<X)r^oDC»o-“C'<n'ta%iot«^cDC»o^c4<n 


B-15 


I 


000 


.866025J7*00  - . 50000007»00  .COOOOOOC  . S66C : 53 '♦00  5000:0 07 -fOO  .00000000 
.9396?260*00  342020 1 9*00  .00000000  . 93969:C'0  *00  -.34202019*00  .00000000 
.99460'74*00  - . 1 7364826*00  .00000000  .96450774*00  -.17364326*00  .00000030 


3 


oooooooooooooooooooooooooooooooooooo 

(O  ID  <D  <0  (D  ID  ir*  ID  (D  (D  <D  <0  (D  (D  ID  (D  ID  (D  >0  (D  (D  iD  (D  (O  CD  (O  CO  CD  (O  (O  <0  CD  CD  (D  ID  CO 

a>  ot  <T>  Ot  (T*  Oi  o>  17'  O)  a»  d*  ni  ct»  O'  f7'  o»  cn  (j>  cr>  o>  oi  di  oi  cti  o»  ni  o>  o*  cii  d C)  oi  o <7> 


O C4fN(Nf>JC>4CNr'4(NOJrN04r'4CNCNr^O4r4r'4fNr4CNCN(N«NCNCNr4CNr4CNr4r4O4CNC4r^ 
t-'  r-  t>.  t-  t-  I'*  t"  t''  f-**  f'  f--  f-  t-  f-*  C'*  C'  I'  f'  I-  I*  I-  t'-  f-»  t-  h*  f-  r-  h-  h* 

c'>fOr'>CDfDrofOrDrDrDfDfnr)">cDro«^cDrofD"7^^fDOror»>ror)ror)r)r)Dror) 

(DcDCDCDCDCDCDCDCDcDCDCOCDCDCDCDCDCDCDCDCDcDCDCDCDCDDCDCDcDCDCDCDCDCOCO 


oooooooooooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooooooooo 

4-4-4-4-f44-  + 4-4-4-444+44-4-44-4-44--f4-4-4-  + 44-+-f4-4-  + 4- 

C>  O)  <7^  Oi  0>  Oi  Oi  O)  Oi  0>  O'  *7>  O'  Oi  (7i  O)  (T)  Oi  0>  O)  O'  0>  CT)  Oi  CTi  Oi  0>  Oi  Oi  D>  U»  Oi  Oi  <7>  <J> 

►1  (N  fN  rv  r4  T'J  IN  f4  f4  rg  CM  04  04  0-4  CN  fN  fN  (N  (N  C4  CN  rv  CN  IN  04  04  (N  (%  CN  fN  (N  O* 

-*  I'  t'*  I-'  t'*  O-  I'  ' t'*  t'-  t - O'*  t-*  t''  C-  I"  I-'  t**  C-  t-*  t'  t'  C'-  I-  f**  t'*  t-«  I"  O'*  h* 

fD  O')  '’I  ro  ^ 'O  O 'T  O i^  o>  o O O O f’)  **)  I’O  f'l  ro  O O)  o>  O O O O o 

ID  CD  ID  ID  'D  ID  <0  CD  CD  'D  iD  'D  'D  ID  CD  (D  cO  'D  CD  iD  'D  <D  ID  CD  iD  iD  ID  CD  D CD  CD  iD  cD  CD  «0  CO 


ID  ID  *D  ID  ID  ID  'O  ID  'D  iD  'D  'D  iD  'D  CD  'D  'D  D (D  'D  -D  CD  D <D  cD  D <D  CD  iD  CD  CD  (D  CD  iD  <D  CD 


CJOOOOOOOOOOOOOOO 
■ I I I I I I I I I I I I I t I 

Oy  04  0 4 r-J  r 4 f4  r4  I"  < 04  C-6  c 4 r>4  f 4 r'4  >"4 

04  -T  '7  ‘T  ” '?  'T  ••  '7  'T  '7  ’T  •?  -J  "T  -• 

<t  X D O)  7/  O;  7 X X X 7 7 Uj  X X D 

•r  'T  '7  '7  -T  '••  -7  'T  '•  'T  ’7  '*  * ^7  *' 

0»  Oi  oi  Oi  Oi  01  O'  oi  Oi  Oi  Oi  Oi  O'  Oi  O'  O' 

t-t** 

(«.  f-  o*  r-  f.  f..  f..  f . I-,  o.  r*,  f*.  f-  o. 


O'TOO'T’T’O'TO’J't'Tnrr^^'T’O’t'J 

oooooooooooooooooooo 

I I I I I I I I I I I I I ) I I I I < I 

C 4 I*  , 0<  CN  ry  r'4  1^4  »'4  04  04  rv  r'4  04  <N  04  04  r'4  04  04 

^7  O -T  'T  *7  O *7  *7  *7  '7  -T  '7  -7  O '7  O O O 't  'f 

o.  o.  Oi  O’  o.  0.  0-  0 o<  Oi  O'  c>  n>  o-  a>  n>  oi  Oi  oi 

X D -71  'J-  T>  '/■  X C X 7v  x.  X r O*  'X  'O  '1.  (D 

-T  ^ 7 '7  '7  -7  '7  •’  *0  O ’7 7 O *7  O 

0>  Oi  Oi  OI  01  Oi  O)  Oi  O'  O)  0»  Oi  Oi  Oi  0»  0>  Oi  OI  Oi  0> 

f-  f'.  r*.  o-  I-'  O'  r-*  r-  f--  o-  f-*  h*  0-  o*  t-  t--  0-  t>*  o* 


f0o>forof0f0f0f0f*)for0f*)f00)f0(oror0(*7<*>o)o>0)rof0f0foroo)r7f0f0 

I I I I t I I t I I I I I I I I I t I I I I < I I I I I I I i I 
OOfJOOOOOCJOOOOOOOOrjOOOOOOOOOOOOOO 
r>roo)roo)r)r)0)0)nnconfnrno)0)i^rorororofOrororofnr)ror)nn 
040404  f-4  04tNry04rN04rNjr404C>4»NO'4r40404f'40404(N04r404C'4CNCNCN04r'4 
04  04f'40404<Nr4  04  04  04  04r'4r'404040404r'404<N04C14(N04r404  04  04  04  04(Nr4 
OOOOOOtDOOOOOOOOOOOOOOOOOOOOOOOOO 


« I i I 
o o o o 
n o CO  o> 

04  CN  C4  C4 
04  CN  04  04 
O O O O 


1000000000000000000000 
I I I I f I ( i I t I ( I 4 I i I I I I I 
’»rr'7*t't'r'T'T’7'T*7'70'Tt0't'T'7'7*T 

>000000  c**oooooooooooooo 

r'-'r'*'t'7'7'7'7'T'7'*'7'^'7'T'7'7'T'T'7»J 
7nr)ncon<or)nr>c^rnr>rnnr)r)r)r)r)0)(*) 
>nnr)r)0)rnroo>r)noir)nfonror>r)fOO)CO 
icDcniniDifiiDcDiDiDiDifiiDioiriDiDininminm 
* I'  I'  I-,  t-  t-  t'  c-  t' t-*  I"  I-*  f'  o*  i'  t'*  o*  t"  O'  o*  r' 
' I'  O'  I'.  I"  I'  O'  t' t'  O'  r-  t-  O'  O'  O'  I'*  j'*  o*  t"  N 


OOOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOOO 
I I I I I I I I I t I I 1 I I I I - I t I I I I I I I I I « I I I I I I 

r4r4  04  04  04  04  04<N04  04  04  04  04CN04r4r4('4r4<N04  04  04  04  04  r40404(N04040404040404 

•D  'D  'D  If  'D  *r  'D  ID  CD  'D  'D  'C  'D  'D  'D  iD  'P  'D  'D  '0  'D  'D  'D  ID  CD  ID  'D  CD  CO  CD  CD  ID  iD  ID  CO  <0 

Cy  f4  f4  f4  f4  Cy  04  04  f4  f4  r4  C4  C4  f4  Oy  Cy  <y  C 4 04  r4  Cy  f4  Cy  04  Cy  r4  04  04  Cy  04  fy  04  04  04  Cy  04 

mnnr'innixxnxfOfO'^xi^fO'^mr'icoxnxif'inxf^rofnnixrinfonfO 

intniDcDiniDinini/ii/iiDincDcDcfii/iiDiniDincntni/iiriniDcfii/iinincniDinminio 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

fonrocoofonf00)nc*)o)f*7ror)00)f0f'>o)o^roo)foof')o)r'7roo>0f'jroo7coo 


’'.ry  04040404040l04CNC404C404C4C4  04  04  04  04  r4040404040404040l04040404040l04 


4 + -f  ♦ 

o o o o 

rwocjo'Docooo 
o o o o •- 

o o o o 

in  m ID  cA 


o o o o 

^-oDoinoo-ooi 
o o o o 

o o o o 

m m in  m 


z •'0<f*>^inioo'CD(Tio^o4r)'vmcro*CD(T»o^o4ro*TincDO'CDOio  — 04r)*TincD 
^•'•-^•-•••-•-•-^(Nfyo4cy(N04  04fyr40<rtr>r>r>r)f*>r> 


••  C4  n V m 


B-16 


00*00009 • 00*00009 


50000*00  .50000*00 


00*00005 • 00*00005 


.50000*00  .50000*00 


♦ 

(/»  ♦ 


S 7 

o o 

Ui  UJ 

O t3 


O ui 

</»  tn 

< D 


O 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o o 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

® O 

m 

•7 

C- 

o 

in 

CN 

*7 

CN  *7 

CN  m o 

»T 

<N 

o 

(L 

f*. 

(/) 

ifl 

'7 

CO 

P) 

PN 

Pi  PN 

Pi  PN  PO 

ro  n 

o 

CO  o 

'T 

”7 

'7 

1/1 

I/I 

Ifi 

C' 

Cj 

1 - 

Ifi 

l/l 

l/l 

'7 

♦7 

n 

•7 

C^i 

CO 

fC| 

CO 

P) 

CO 

CO 

ro 

P)  CO 

P) 

ro 

ro 

ro  ® 

»- 

o 

O 

o 

O 

O 

CJ 

o 

O 

o 

»- 

o 

o 

CJ 

o 

o 

CJ 

Cj 

Cj 

O 

o 

Cj 

o 

Cj 

<j 

O 

fj 

O 

o 

o 

o 

o 

O 

O 

o 

o 

o 

O 

o o 

o 

O 

o 

O O 

o 

C'J 

CJ 

o 

CJ 

O 

O 

O 

CJ 

CJ 

Cj 

CJ 

Cj 

o 

O 

o 

o 

o 

o 

o 

O 

O 

"• 

o 

o 

o 

01 

Ot 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

O 

O 

o 

CJ 

o 

o 

o 

o 

o 

o 

o 

o 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

O CN 

'T 

»>. 

O 

in  (N 

*7 

CN 

»7 

CN 

in 

o 

t>- 

»r 

CN 

O 

f- 

in 

in 

cT 

P) 

ro 

PN 

Pi 

PN  Pi 

Pi 

ro 

P> 

*7 

in 

OJ 

n 

o 

»7 

o 

If- 

i/) 

10 

1* 

Cj 

{ • 

«fi 

in 

m 

•7 

'7 

»7 

'7 

CO 

fO 

rc; 

CO 

CO 

Pi 

CO 

ro 

Pj 

r'l 

ro 

ro 

ro 

CO 

fO 

p> 

CJ 

o 

O 

o 

o 

CJ 

O 

CJ 

CJ 

cj 

ij 

Cj 

O 

CJ 

r > 

CJ 

CJ 

CJ 

o 

o 

Cj 

CJ 

Cj 

o 

O 

Cj 

o 

o 

CJ 

Pj 

PJ 

o 

o 

o 

fj 

o 

CJ 

C) 

fJ 

rj 

o 

CJ 

o 

CJ 

CJ 

O 

O 

rj 

O 

o 

O 

o 

o 

O 

CJ 

Cj 

PJ 

Pj 

o 

Cj 

Cj 

o 

o 

O 

o 

"■ 

o 

o 

O 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

O 

O 

o 

o 

O 

o 

o 

O 

O 

O 

CJ 

O 

O 

O 

O O 

o 

o o 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

■ 4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

* 

4 

4 

4 

4 

4 

4 

4 

4 

C'i 

n 

t • 

o 

in 

CN 

•7 

CN 

'7 

C4  iTi 

o 

r~ 

"7 

CN 

O 

f * 

Ifj 

ifl 

'7 

ro 

fO 

Pi 

Pi 

Pi 

Pi 

Pi 

Pj 

CO 

•7 

m 

m 

•T 

*T 

•7 

if- 

ifi 

'0 

t - 

fj 

I" 

•£» 

ifi 

/- 

'7 

•J 

'T 

•7 

ro 

CO 

fCj 

P| 

ro 

ro 

P| 

CO 

CO 

ro 

ro 

f'. 

fO 

ro 

(0 

fO 

CO 

O 

O 

O 

CJ 

Cj 

o 

o 

«- 

CJ 

tj 

Cj 

( J 

Cj 

Cj 

Cj 

Cj 

o o 

C'J 

O 

o 

Cj 

Pj 

c J 

CJ 

O 

O 

c J 

o 

Pj 

CJ 

o 

PJ 

o 

o 

O 

o 

o 

o 

o 

CJ 

CJ 

O 

CJ 

CJ 

CJ 

O 

CJ 

O 

o 

O 

O 

Cj 

CJ 

CJ 

O 

CJ 

O 

CJ 

rj 

O 

PJ 

PJ 

CJ 

PJ 

o 

"■ 

o 

o 

rj 

CJ 

r j 

o 

C J 

rj 

»J 

r> 

CJ 

CJ 

o 

C J 

O 

CJ 

c> 

f J 

CJ 

o 

o 

Cj 

r J 

o 

CJ 

PJ 

PJ 

PJ 

o 

CJ 

PJ 

o 

PJ 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

•T 

t* 

O 

I/I 

»7 

C4 

•7 

C4 

l/l 

CJ 

1 - 

•■4 

CJ 

cf) 

1 

l/l 

l/l 

'7 

CO 

f, 

Ci 

r i 

ri 

CN 

Pi 

ro 

CO 

'7 

If) 

m 

!>• 

'T 

'7 

1/. 

i/i 

■n 

1 - 

CJ 

1 - 

>0 

i/i 

in 

•7 

'7 

•T 

•1 

f*’J 

CO 

CO 

ro 

CO 

CO 

f'i 

•'i 

ro 

ro 

CO 

10 

CO 

CO 

P> 

CO 

ro 

(O 

O 

f j 

CJ 

CJ 

CJ 

o 

•- 

t J 

CJ 

1 J 

*j 

CJ 

CJ 

O 

Cj 

CJ 

o 

». , 

Cj 

O 

r J 

fJ 

< J 

J 

o 

( J 

rj 

C J 

Pj 

CJ 

CJ 

C 1 

o 

O 

O 

O 

O 

CJ 

CJ 

o 

O 

o 

CJ 

o 

O 

CJ 

O 

f J 

CJ 

r J 

o 

'.J 

o 

Cj 

fJ 

CJ 

CJ 

( J 

CJ 

CJ 

Cj 

CJ 

PJ 

CJ 

rj 

o 

CJ 

O 

o 

O 

o 

Cy 

Cj 

C/ 

CJ 

CJ 

CJ 

Cj 

o 

O 

o 

Cj 

CJ 

O 

CJ 

Cj 

f . 

r J 

C J 

o 

t i 

PJ 

Pj 

r • 

CJ 

Ci 

O 

o 

* 

4 

t 

4 

4 

4 

• 

4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

* 

» 

4 

4 

* 

« 

4 

4 

4 

4 

* 

f-* 

o 

I/I 

r 4 

•7 

r 4 

’7 

C 4 

If. 

C J 

t • 

n 

r^ 

C J 

1 - 

if. 

in 

•T 

'O 

Pi 

r 4 

C . 

r 4 

f , 

CO 

CO 

• 

lO 

Ifi 

1 * 

m 

»r/ 

'7 

1// 

1/. 

1 • 

Cj 

f • 

0 

.f- 

If. 

’J 

•7 

• * 

'» 

''j 

cc- 

'0 

CO 

'0 

"j 

"I 

"■< 

p< 

•'t 

rO 

CO 

CO 

ro 

Ci 

o 

CJ 

c> 

Cj 

•- 

» V 

• J 

CJ 

Cj 

CJ 

o 

- ^ 

o 

<J 

1 J 

< t 

C - 

o 

( J 

• J 

o 

« J 

' J 

Cj 

•J 

1 • 

f > 

' J 

o 

o 

o 

o 

O 

CJ 

rj 

O 

CJ 

CJ 

Cj 

rj 

O 

o 

CJ 

Cj 

CJ 

o 

<j 

O 

o 

Cj 

' . 

rj 

Cj 

< i 

<j 

c y 

fJ 

C J 

fj 

o 

♦ 

f 

4- 

♦ 

f - 

•- 

•- 

*■ 

*" 

•" 

•" 

•“ 

- 

- 

*" 

*■ 

•" 

•- 

*" 

7.  O 

■f 

UJ 

o 

o 

♦ 

(U 

u. 

♦ 

<D 

i< 

•f 

o o 

o 

o o 

O 

tj 

o 

o 

O 

O 

o 

o 

O 

O 

O 

o 

o 

o 

o 

o 

Cj 

O 

o 

O 

o 

O 

O 

o 

O 

O 

O 

o 

(T 

(/) 

4- 

4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

♦ 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

< 

o 

+ 

o 

in 

r4 

•7 

C4 

CN 

in 

O 

r • 

'T 

CN 

o 

in 

f» 

m 

in 

'7 

CO 

ro 

f'i 

CN 

f i 

Pi 

Ci 

CO 

ro 

•7 

in 

in 

r> 

® 

O 

(r 

♦ 

ir» 

m 

•n 

I- 

O 

(' 

O 

ir> 

l/l 

•7 

»7 

'7 

'7 

CO 

fO 

fO 

ro 

C'l 

CO 

CO 

'O 

fO 

P| 

"i 

Pj 

p) 

pj 

Pi 

CO 

ro 

(O 

to 

•7 

O O 

f) 

4 

o 

o 

o 

CJ 

r' 

O 

'J 

Cj 

o 

CJ 

Cj 

O 

CJ 

CJ 

o 

<j 

CJ 

O 

o 

CJ 

Cj 

< J 

C J 

O 

CJ 

«J 

o 

<j 

C J 

(J 

CJ 

rj 

O 

< 

4 

o 

o 

o 

o 

CJ 

O 

O 

CJ 

o 

O 

CJ 

O 

O ?J 

CJ 

CJ 

o 

CJ 

CJ 

rj 

CJ 

CJ 

cj 

o 

CJ 

o 

Cj 

PJ 

c> 

CJ 

o 

O 

O 

UJ 

4 

O' 

o 

4 

• 

• 

o 

4 

® 

4 

(I) 

4 

4 

o 

o 

o 

o 

rj 

o 

CJ 

CJ 

o 

o 

O 

O 

a 

CJ 

o 

o 

o 

o 

o 

O 

CJ 

CJ 

CJ 

o 

O 

CJ 

O 

O 

o 

o 

PJ 

CJ 

o 

♦- 

® 

4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

♦ 

♦ 

CO 

4 

in 

C^ 

o 

CN 

♦7 

CN 

in 

CJ 

r- 

•7 

CN 

o 

Oi 

1 ' 

in 

m 

•7 

CO 

CO 

Pi 

Ci 

f«l 

Pi 

Ci 

ro 

ro 

l/> 

m 

I'- 

ro 

O 

PN 

S 

in 

4 

in 

•I> 

t - 

O 

t - 

'0 

I/I 

1/1 

'7 

»♦ 

o 

'^1 

to 

CO 

CO 

CO 

C'l 

r^l 

CO 

CO 

C'J 

'0 

Pl 

CO 

CO 

P7 

CO 

CO 

ro 

ro 

'7 

't 

a) 

4 

o 

o 

CJ 

>- 

O 

rj 

o 

CJ 

o 

CJ 

CJ 

o 

o 

CJ 

o 

o 

CJ 

Cj 

CJ 

o 

o 

rj 

CJ 

PJ 

CJ 

PJ 

PJ 

PJ 

o 

o 

O 

CJ 

CJ 

>x> 

• 

4 

o 

o 

O 

o 

O 

O 

o 

O 

o 

o 

CJ 

CJ 

o 

o 

CJ 

o 

o 

Cj 

o 

CJ 

o 

cj 

o 

PJ 

CJ 

O 

O 

O 

o 

o 

o 

O 

CJ 

m 

V- 

»- 

»- 

r* 

r* 

w- 

w 

•— 

o 

4 

• 

• 

►- 

r** 

_i 

4 

te 

CO 

3 UJ 

4 

UJ 

n 

u. 

re 

4 

7 

4 

•- 

•- 

•- 

•- 

— 

•- 

w 

•- 

w- 

•- 

•- 

— 

*- 

»* 

— 

•- 

•- 

4- 

•- 

4- 

o 

o 

4 

in 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

CJ 

o 

o 

o 

o 

o 

o 

O 

o 

o 

CJ 

o 

o 

o 

o 

o 

o 

o 

O 

4 

Ui 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

t- 

<N 

o 

CN 

'7 

CN 

in 

rj 

f - 

»7 

CN 

CJ 

(D 

f* 

m 

in 

•7 

CO 

CO 

Pi 

PV 

CN 

PV 

PN 

CO 

P7 

•7 

in  m 

/- 

© 

O 

Pi 

't 

u r> 

4 

«t 

CO 

t- 

O 

I*. 

CO 

I/I 

in 

n 

•7 

^7 

f*t 

CO 

CO 

CO 

CO 

CO 

fO 

CO 

CO 

CO 

ro 

ro 

P| 

P> 

P) 

p) 

CO 

CO 

ro 

'7 

*7 

*7 

n n 

l/> 

4 

7 

ij 

CJ 

•> 

O 

o 

rj 

CJ 

CJ 

CJ 

Cj 

CJ 

CJ 

o 

Cj 

CJ 

CJ 

CJ 

CJ 

Cj 

o 

O 

Pj 

CJ 

o 

CJ 

O 

o 

O 

o 

o 

PJ 

O 

O 

(N 

Ui 

4 

M 

o 

o 

o 

o 

o 

O 

CJ 

O 

o 

0 

0 

CJ 

o 

Cj 

o 

o 

O 

CJ 

O 

o 

o 

O 

o 

o 

O 

O 

o 

o o 

o 

O 

o 

O 

o 

(j 

4 

a 

•- 

•- 

»- 

•- 

•- 

•- 

•- 

r- 

»* 

»- 

w- 

•- 

•- 

»- 

♦- 

•- 

4- 

4- 

4- 

u 

4 

a 

• 

• 

• 

> 

• 

4 

o 

UJ  i/i 

4 

o 

o ® 

O UJ 

4 

o 

Ui 

(N 

4 

V* 

»- 

•• 

w- 

«— 

w 

4“ 

4- 

4- 

(/) 

UJ 

4 

? o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

0 

c 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o o 

cr 

o « 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

nr 

« 

4 

UJ  D 

n 

CN 

n 

r4 

m 

D 

•> 

'N 

n 

f*^ 

t' 

in 

IP 

*^ 

CO 

CO 

CN 

Pi 

CN 

Ci 

PN 

ro 

ct 

m 

m 

P* 

® 

O 

PN 

^7 

r>4 

a 

a 

4 

O -t 

C4 

1 

O 

i • 

'Ii 

ifi 

I/I 

•7 

"7 

'7 

•7 

rc| 

CO 

C', 

n-i 

CO 

CO 

CO 

ro 

C'l 

P> 

ro 

C'l 

CO 

ro 

ro 

p) 

ro 

ro 

P) 

•^7 

»7 

•7 

c 

O 

3 

O *t 

4 

• u 

o 

Cj 

o 

CJ 

Cj 

CJ 

O 

O 

O 

o 

o 

o 

O 

Cj 

CJ 

o 

o 

Cj 

O 

CJ 

o 

CJ 

O 

O 

o 

o 

o 

o 

O 

O 

O 

o 

< 

4 

o 

o 

O 

o 

o 

O 

O 

O 

O 

o 

o 

O CJ 

O 

o 

o 

CJ 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

O 

O 

PJ 

U ® 

4 

> Z 

•- 

•- 

•- 

•* 

w- 

»- 

•- 

w- 

•- 

•- 

w- 

•- 

»- 

•- 

Ui 

CO 

4 

u — 

• 

* 

* 

• 

o< 

4 

X 

4 

D 

4 

a 

•- 

»• 

♦- 

»- 

•- 

•- 

44 

• 

4 

►- 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

0 

0 

0 

o o 

o 

o o 

UJ 

y 

tn 

4 

•f 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

* 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

o 

a 

4 

jt 

(N 

o 

r^ 

m 

r> 

f- 

o 

Ci 

C) 

® 

» 

if. 

u. 

•7 

CO 

CO 

CM 

f 4 

CN 

f'i 

Pi 

CO 

CO 

♦7 

l/l 

l/l 

® 

o 

PN 

<7 

O 

Ui 

O 

4 

O 

C-* 

'i‘ 

u'l 

*7 

•7 

o 

<*•1 

ro 

n 

'O 

CO 

CO 

'*1 

CO 

ro 

CO 

CO 

CO 

CO 

CO 

ro 

CO 

p) 

P)  CO 

•7 

•7 

«7 

*7 

m 

;? 

O 

4 

m 

». 

o 

o 

O 

o 

O 

CJ 

O 

o 

o 

O 

o 

o 

CJ 

o 

CJ 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

O 

O 

O 

o 

t 

• 

Cl 

4 

•ft 

O 

o 

o 

Cj 

o 

O 

o 

o 

o 

o 

o 

Cj 

o 

CJ 

o 

CJ 

o 

o 

o 

o 

o 

o 

o 

Pj 

o 

O 

o 

o 

O 

o 

O 

O 

o 

►- 

4 

UJ  < 

«• 

•* 

•- 

•- 

•» 

4- 

4- 

*< 

► 

4 

UJ  s 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

♦ 

• 

O 

O 


♦ ♦♦♦♦♦♦ 


B-18 


k. 


ooo  oooooooooooooooooooooooooooooooooooo  ooooooooooooooooo 

♦ + + + 4444444444444444^ 

if>  i/»  I • CN  < 4 <n  «7  iTj  i/i  t • a>  o f4  *7  J ' o i/i  C4  *7  r4  *r  tN  i/>  o I-  n r4  <J  CD  f"  i/»  ifi  *7  po  r4  cn  o cd  r-  in  m '7  o ro  cn  (N  f4  C4  fO  fo  ^ 
n o ro  ")  po  rn  ^ m n n po  *7  *7  n *7  m u>  <ly  i • o i - »d  i/i  i/j  *7  '7  *7  *7  r)  rj  cn  f*)  c*)  o *7  '7  rn  n ro  c'j  r')  n n m r>  f*)  r)  <*)  r)  n 

OOCjOtyC/OOOorjOOOOOOOU^iOOO’-fJOOOO  JOOOOCjOOOOCjfjOOOOCiOOOOOOOOOOOO 

ooo  ocjooooooooooooooooooooooooooooo*;^oooo  oopjoooooooooooooo 


ooo  oooooooooooooooooooooooooooc>oooooooo  ooooooooooooooooo 

44+  ♦+444444  + 444444444444  + + 4444  + 4 + 4»  + 444  444  + + 44  + 4 + + 444  + 44 

mt->no  C4r4(NforO'7ir>ir)h'CDor4<7t^oinr4»7r4'7r4mor-«7f'4oajr~inm*7fnnr4f4  ocDt''mi/>'7nr0fMC4f4O4O4(*)(*)«7i/> 

n n ro  'll  rn  m f*)  O '7  *7  *7  *7  t/»  a.  tl>  I • o f'  ID  iT*  if»  *7  •’7  *7  •)  m r»)  ro  cn  O c')  cn  C7i  «7  ro  fO  c*>  fo  ro  ro  co  n ro  ro  cO  fO  c*)  O ro 

O O Cj  •-  O O <>  O O O O O O O O O O O O O <-»  •-  CJ  O O O O O c>  CJ  O <J  O O O 1-#  O O C4  O O O O O O O Cj  o o o o o o o 

ooo  o o <■>  (i  o o < j o o o o o o o o o o o f7  C7  o o o f > fj  o r j f-5  o f > c j o o o rj  o ooooooooooooooooo 


ooo  ooooooo  ooo  fDoooooooo  ooooooooooooooooo  OrjOOOOOOOOOOOOOOO 

444  444»4444444+444444444444444444444444  4444444444444444  + 

r*-  o)  o tN  n fo  *7  iT)  in  f-  OD  o r4  *7  f • o If)  74  *7  <N  *7  i/i  o I-  ^7  oj  o CD  f-  i/i  ifi  *7  r4  04  CD  f-'  tr*  in  a o)  cn  74  04  CN  04  74  <n  f*>  *7  in  in 

70  T'i  ■?  CD  o,  f*,  fO  frj  fr,  (r^  (»-,  r*  .y  j,  tl*  <£.  I - 7j  f ' '£>  l/i  if/  *7  •‘T  '7  '7  70  fO  7^  70  (')  r»-,  rr,  f<-,  rr)  r-/  Qj  f»)  /▼,  f»/  f*!)  fO  T'/  70  70  70  70  ro  fO  T')  7*>  70 

7j  <>  O •-  O O O <>  <J  <-»  7_,  Cj  o 7>  tv  O C«  O 7J  O ^ O O O 7j  o O O O 7^7^00  fj  Oj  Cv  0 0 7,0  tv  o O 7J  ♦>  O 7j  7>  rj  o O c.  C.J  o O 7J 

O 7 J O 7J  O fj  O O O O 7J  C>  O 7>  o O O O r>  o O O O O O O O 7j  o O O tj  O O O O O tj  O O O O O 7J  O O O 7J  O O 75  o OJ  o o o 


O O C5  O C>  7.5  75  t5  O O O O « 


7>  75  7v  o O O 7j  <j  o O 05  7>  75  o O O O O O 75  7_>  tj  o . O O C5  O C>  7.5  05  t5  O O O O « J 05  O 05  75  05  O 05  75  75  O 75  75  fJ  O O5  O 7j  O 
444  44*4*44t«4444444444l  .444444444444’4«4  444444444444444  + + 

O)  75  74  70  70  *7  ifl  iTM  ' OD  O 74  'T  » ' O iTl  74  *7  74  CT  74  in  fj  f - *7  7-4  o to  | - 1/,  ifi  • 7 70  'O  74  74  7 4 74  f-  l/i  lf»  7»j  7»)  04  O4  74  74  74  70  70  *7  ifi  ifi  h' 

O -7  *7  I ' 70  r^>  70  70  'O  70  70  'O  *7  '7  '7  '7  l/»  U»  'Ij  » - 7J)  l • O i/i  lOi  ’7  '7  *7  'I  <0  70  70  70  70  70  70  70  <0  'O  70  I 70  «0  'O  70  to  to  70  70  to  T'j  7»)  (O  70  70  70  7*) 

7./  ^ 75  7_/  f , <->  O 75  t5  7j  <5  •'V  7j  (5  I 5 O O t5  »-  O t5  O C5  O O C5  O » 5 O < 5 O 75  7.5  t5  ( 5 ' j T.j  C5  7 4 < 5 75  < V O ' 5 O O Cv  O ' > ( 5 7J  r5  O O O O 

75  O ‘5  0 0 0 75  t5  00000  75  75000000  7j0000000f5  000000000  O 75  O 75  0000  75  ooo  00000 


> O 7-5  75  7^  75  O 7/  t>  7 5 rj  « 


7 > O O O O O O O O O 75  o t 5 75  O O 75  O C>  7v 


O O O O O O O O 75  O O 0 O O O O O 


/ 74  *T  70  70  * 7 . 71  ifi  r - X 7-5  74  T « - O m 7,  'T  74  *7  74  UJ  O f ' «T  74  rjcfjt-  i/»  iT.  '7  t^/  to  r.,  7 4 7 j 7 4 7 4 lT.  i/,  *7  7<,  *0  74  74  74  74  74  ' j ^o  *7  i/<  iT*  f"  X 

•;  *7  T 'D  f'i  T')  /»,  ro  'O  "j  'r  ‘7  'T  i7<  .7/  ■£>  I ■ O t - X 7«  J»  ^7  *7  '7  <0  to  f ro  fO  X f',  70  'o  fO  '0  'O  *0  f'<  70  X t*-,  'o  fO  'O  X 7*4  7n  to  70  fO  'O  70  70 

'5  «■/  <5  »-  75  l_5  C>  ' y fj  '.V  0 7-5  O ‘ / O 75  fv  7 »»-  O O O 75  O O <5  «->  '5  o 75  O 7>  .>  «5  «_/  f >»/  7,  Cv  O O O .->«  5 ■ ‘-V  <5  «-)  O Cv  C5  '5  7 > O O 

75  00  0000  75  75  0000000  7j0<j7jo000t50t/0t5000  75  7.5  c5c50o0  0«-50<-y00t5  0 75  000000  75  0 


OOO  ooooooooooooocoooooooooooooooooooooo  ooooooooooooooooo 

♦ ++  +♦♦♦♦♦  + ♦ + + + ♦ + + ♦ + ♦♦  + + + + ♦ + +♦  + ♦+  ♦+  ♦♦  + + + ♦♦  + ♦ + + ♦ + + ++  + ♦+  + + + 
f j «7  r-  70  '7  in  in  f ' X o 74  ^ c-  o in  74  r4  77  r4  in  75  I"  77  04  o CD  t->  i/i  in  »7  70  70  C4  f4  04  74  r4  70  m n 70  70  74  74  74  74  04  70  70  '7  in  in  o>  CO  o 

T *7  *7  in  7'’.  70  f'i  tt,  r*)  P*  '7  -7  ••7  -7  1/,  X f Cj  t-  •!>  l/i  if»  '7  '7  '7  '7  fO  70  fO  70  70  70  <0  7<)  7O  70  7Or74f»>i/,7n<Of'J'0  7O70f0  77/P)7O7O7OrO(0r0  70»T' 
75Ot>^75OO75OC5O75(5Ot->'j75O»-05OOOOOOOt5OOOOOOOOOOOO(4OOOOO  75  75  C5OOOOOOOOO 

75  00  7500  0 00  0000  75  00000000000500000000  0 00000  OOO  t5  000000000  0000 


OOt5  OtJ  75  C5  05  00  00000  05  05  0000000000  00000  000000  05  C/OOOOOOOOOOOOOOOO 

♦ ♦♦  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦4  + 4 + + + * + + + 4 + 4+  + + * + + *+  ♦■  + + 4 + 4 + + + 4 + + + + + + + 

•t  f-  o *7  I/)  in  C'  X 75  74  *7  r-  o I/I  74  *7  <4  *7  04  ifi  75  I-  '7  04  o no  I'  in  in  '7  7n  7*7  74  7N  74  74  04  <7)  *7  7*)  70  74  74  74  74  04  x 7n  *7  in  m x o m 

*7  'I  I/I  ^7  70  'O  »'»  'O  X *7  '7  ’7  '7  i/i  i/i  </>  t'  C5  I - <0  ifi  l/i  ‘7  *7  '7  '7  X X 70  fO  7*)  70  70  X X X X X 70  X *7  X 70  X 70  X fO  fO  X X X X »0  X X X <7  <7 

OO<5»-OOO75rjoOOOO75O7L»*-O05Ot5f5OOOOOr5O75r5  75  0>OOOO7jor4OOC5OOOOOOOOOOOOOO 
Cj  rj  n OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO  00O000O0O  75  7JO00000 


ooo  oooooooooooooooooooooooooooooooooooo  ooooooooooooooooo 

♦ ♦♦  ♦+  + ++  4 + + 4 + + + + + 4 + 4 + + + + + + 4 + + + 1 + + + + + 44+  4 + 4444f4+  + +4  + + + + + 

f-75in  ininr'*xox'7f  0i/»f4«774»7<‘4i/i0r'*«4r40xe'-inin»7xxr4r4xr4  74xx*7  xxr4r4r4r4C4XX^ininh*xor4«T 

•7  i/i  I/I  X 70  70  X 70  X '7  '7  '7  1/1  ifi  X I - O f - X I/I  l/i  *7  '7  *7  '7  X X X X X X X X X X 70  X X X 70  X X X X X X X X X X X X X X X »7  *7  ^ 
0 75  0*'Of5  75  0 00  0 0 7^<5  0T5-OOt50r5  000r5  00f5  00  0 75  0000000r4000  75  0000000000000 

ooo  oooooooooooooooooooooooooooooooooooo  ooooooooooooooooo 


ooo  oooooooooooooooooooooooooooooooooooo  ooooooooooooooooo 

♦ 44  44444444  + 44444444444  + 44  + + 44  + + 4+  + + + 44  + + 4 + + 44  + + + + + + + + + + 

Oi/lTN  in»'CD7~>f4*7r  0inP‘*7X'7XinO7'*-'7XOXf^inifi'TXXr4r4r4r4r4XX'7in  Xf4f47Nr474XX»7inin+^CD07N«7r^ 

m l/i  ID  f4  X X X '7  n ^7  n 1/|  iTi  ID  7'  O I i3  i/i  ifi  ^7  *7  '7  '7  X X 70  X X X X X X X X X X X X X 04  X X X X X X X X X X X X X »7  «7  *7  «7 

o 75  o»-o  7_,  75  oo  ovooooc5»-  7-5  00000  7.5  ooo  Ofjo  75  oooooooooo  74  ooooooooooooooooo 

ooo  oooooooooooooooooooooooooooooooooooo  ooooooooooooooooo 


ooo  oooooooooooooooooooooooooooooooooooo  ooooooooooooooooo 

4 + 4 ♦444444444*444444*444444444t44  + 4 + *44  +4*4444+4++4  + + + + + 

<n  74  *7  r-  CD  o 74  »r  t - o m r4  *7  74  n '■4  iTi  o I *7  74  o cd  r-  i/i  i/i  n x x 74  74  tn  74  *4  7»»  x *7  in  i//  74  <4  74  74  x x x *7  in  in  co  o <4  *7  o 

in»Dt'  ^ XX  n X X ift  tfi  <0  t-  cj  t ■ 'Dini/ix  *7  x x "jxxxxxxxxxxxxxxxx^-xxxxxxxxxxxxx  xx^in 

0 0 0«“0000000000*-  000000  0000  O OOOOOOOOOOOO  O.  07400000000000000000 

ooo  00000000000000007500000 oooooooooooooo  ooooooooooooooooo 


• <Nnv^n^£)^•eo<^o-r«lX4Tln^o^•®<J>o  — rnx^mipr^ooo'o  — r^x^Tiftco 

•-••»'^*-»-*-^»-*-f4Mf4Ci<<4r#f4<4f4C4XXXXXXX 


ooooooooooooooooooo 

inif>r^  aDocN«ir>ou^f^^  <N«TfMU^ot^^ 
ro-*)'nfO»T*7n»tifiir»«£tt^or'ir)inin^^ 
ooooooooo  ooo  --oooooo 
ooooooooooooooooooo 


ooooooooooooooooooo 

ini~a)ofN«T  h-otfJCN«tfN*i<Ninor'’»ci 
'T'7 

00000000000--0  000000 
ooooooooooooooooooo 


ooooooooooooooooooo 

t'CCiOfN'Tf''Oir)CN'TfN*T  CNU10h»«J<NO 
r*i  n »I  t 'T  1 iTi  iTj  'O  t--  O f**  'ii  in  *7  ^ *T 

oooooooooo  — OOOOOOOO 
OOOOOOOOOOOOOOOOOOO 


OOOOOOOOOOOOOOOOOOO 
* * * 

cDCir^«7i-oir»r^»T<N'7r<inof-'*TfNO0D 
f'J  n *1  'T  'T  U»  I/*  'i;  f • O I • 'O  i/<  i/i  'T  'T  'T  't  ro 

O O Ci  rj  o O CJ  O O ^ O O O O «J  o o o o 

OOOOOOOOOOOOOOCiOOOO 


O O O o O O O O Cv  o O O o f)  O O O O O O O O O O O O O O O O O O O O o o o o o o o o o o o o o o o o o o o 

O r4  'J  (-  O I/)  *7  <■-*  n i/i  O t-  *7  O 'n  f-  *-1  C4  \S,  O f-  *7  C*  O fTj  I'  Ji  il>  ’1  ro  fO  r-^  O Cv  rv  «■<  »7  iTi  ifi  t"  CD  O '7  r • O i/)  ra  «7 

'T  *7  *7  if  It  >0  i'  r>  f'  'C,  l/i  l'.  •?  '?  •!  '7  m 'i  t - ‘i>  ■.{,  Jt  'T  n -7  ('I  ")  ")  r)  fn  'n  ">  ro  ^ *7  -7  n *7  i/J  in  <n  f*-  o 

O O O O O — O O 1_»  O O O O O O O f)  O O O O CJ  O O Cj  '■J  cj  O O O O O O «■«  OOcjoOOoOO<’jOOOOO<’><'>»- 

O O O O O O O O O O O O O O O CJ  O O O O O O O (_»  O O O O O C>  O O O O O O O O O O O O O O O O O O o o o o o o 


OOOOOOOOOOOOOOOOOOO  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

44444444444444444-f4  +-4  + + -f444444444-f444444  + + 444-4444444'f4--f 

ry»7f-oin<N*T(N'7r^inof'-'7<NOcor^in  oinot''Tf^ooot'-inin'TrofDfNfNfNfN<NfnfO*74ni/»t^aoo(N^r>‘Oinr>i‘tr<«7 

'T  '7  *7  i/.  i/.  '£»  f-  O J ' 'X>  if*  U»  t *7  *7  '7  O ro  in  'D  iTj  i/J  't  ■‘7  'T  '7  CD  O fO  ro  ro  fO  o 'n  O '■n  f')  ro  o o m ro  O ^ *7  'T  *7  iTi  in  '£*  f'*  o r* 

oofJoooo'-ooooooocjooorDooooooooooooooooooooooooooocioooooo»-o 
ooooooooooooooooooo  oooooooooooooooooooooooooooooooooooo 


ooooooc:>oooooooooooo  oooooooooooooooooooooooooooooooooooo 

♦ 44444444*4444444f4  4444444444444444444444444*4444444-f44 

•Tt'Oif»r*«7<N'7c-4inof''7<Noa)j'*i/'tin  inor'‘'7(NocDh*inin'7f*icr)c>jrNr’^r4fNrocn«7ii“iinr^a5  0f’<«7i'-oirifM'7CN«7r^ 

'T  '7  l7<  i/*  '0  » * O f - 'D  l/»  t/»  '7  ^7  *7  'T  CO  m CD  O *7  i/>  in  't  *7  "7  T CD  CO  fn  CO  CO  CD  CD  O CO  c-^  r>  CD  n CD  'D  cn  CD  cn  'T  ^7  *7  T If*  in  'fi  I'  O r'  CD 
OOOOOO  — 0 000cj000<j00  0c»>0  00  00000000000000f->000000000000000»-00 

ooooooooooooooooooo  oooooooooooooooooooooooooooooooooooo 


ooooooooooooooooooo  oooooooooooooooooooooooooooooooooooo 

4444444444444444444  ♦♦44444444-44444--f4  + 44+444444444'f444'f4- 

r-ojnr^'7r^'7r4ino4^»7tNOcor^tnin'7  of^*7<Noa)h-ir)in'7c*)CDrN(NfNc>4fNCDCD*7ininr-'a3oc>4'7f'»oincN»7f><^cNin 

•7if»if»»D»'OI''Dif*in*7'7'7'7CDfDCDCOfnf*)in'7'7't*7c*)CO(DCDf*)nr)rDfDfDCDr)rDCDrDrDf*)f*)Ot*7*7'7ininiD«’*Or^(Dm 

ooooo--oc:>oooooooooooc»^»uocjooooooooooooooooooooooooooooo»-ooo 
OOOOOOOOOOOOOOOOOOO  oooooooooooooooooooooooooooooooooooo 


ooooooooooooooooooo  oooooooooooooooooooooooooooooooooooo 

♦ ♦♦■♦•f444-44444  + -f4-444  44444444444444444444444444444^  + 44-44-4  + 

oin<NT<N*7r4ir>or‘*tfNOT5r'.inin*7CD  f-'7rMOcnr-inin'7CDr)r4rMC<<Nr'<cDCD'Tinini-a)OfN'tf-oino*<7f'4'7«Ntno 
inir»‘Dt-Of  *Di/tin'7'7'T'7CDcor>OCDr><N«t'7'7'Tc«7cDr)fDCDrDOr)r>f*>CDCDCDCDfDf*)CDfnf*)'7'7'7*7inincOf^Or'tDinin 
0 0 00*-0000000000000oc*)000000000000000000  0 0 00000000000*-0000 
ooooooooooooooooooo  oooooooooooooooooooooooooooooooooooo 


ooooooooooooooooooo  oooooooooooooooooooooooooooooooooooo 

44444444444+4444444  444444444444444+4  + 444*444444444444-44 

inr4*7r4’7<Ninoh-'7nocof'ir»ir)'7f*>c*>  »tr4oaD7'inin«7'*7cn<Nr«fStN(Nc«oco»7ininh‘a>or4»7h*onc'(*7<N«7o<»nof^ 

ln'D^•0^''Ol/■»ln'7't•7■'7cncDCDHc^c*)^»)^'J»7'7^Dr»)fr)^^fr)fr)rl->r+)r»>f<»)^»^r)r»)'DCDrDf*)r>co•7'7'7•7lnncOf•-Ot'•'Dlnln^ 
ooo*«ooooooooooooooocDoooooooooooooooooooooooooooooo«-ooooo 
ooooooooooooooooooo  oooooooooooooooooooooooooooooooooooo 


(N  r)  2 <0 

o f*) 


f 


o 

o o 

»- 

o 

O 

o 

o 

o o 

O 

o 

o 

4- 

o o 

o o 

— o o 

o 

0 •- 

o 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

Cj 

o o 

o o 

o 

o 

o 

o 

o 

o o 

o o o 

0 0 

•f 

1 

1 

) 

4 

X 

X 

X 

X 

X 

4 

4 

X 

X 

4 

4 

4 

4 

X 

X 

X 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

4 

X 

4 

X X 

X 

X 

X 

X ♦ 

o> 

o> 

o 

ro 

(J>  cx 

n 

L'l 

'D 

ID 

O 

o 

o 

X r- 

«- 

in  ID  rx 

«- 

o 

o 

t" 

X 

If* 

'D 

X 

44 

ro  cx 

!■«. 

X r* 

ro  X (D 

<0  X 

o 

o 

o 

h' 

o 

rr- 

o 

n 

;fl 

o 

CO 

ro 

ro 

9' 

'D  ID 

9 

I- 

ID 

•7 

f- 

o 

If) 

If. 

w 

X 

CN 

X 

X 

4- 

m 

ID  9. 

X 

f'- 

■D  ID 

ID  O X 

o CD  in 

•- 

<X)  rx 

rx 

n 

® 

•- 

0. 

'0 

•- 

o 

0. 

0i 

rx 

O 

3' 

9 

o 'D 

o 

ro 

in 

9 

I'- 

9) 

X 

4- 

o 

9. 

9. 

rx 

4- 

0* 

9 

O ID 

O 

X - 

ro 

9 

C'* 

cx 

CD  9 

O. 

(£> 

CTj 

O 

CO 

CJ 

ro 

rx 

•- 

If. 

rx 

ro 

o 

•- 

cx 

rx 

(T) 

9 

ro 

f-- 

X 

o 

CO 

o 

X 

cx 

4- 

If.  rx 

X 

O 

4- 

fX 

rx  CD 

9 

X n* 

X X 

'’• 

(y» 

ID 

U'l 

1 

,* 

1 

,* 

9. 

CO  ro 

1 

1 

4 

CO 

i 

f ' 

9. 

X 

1 

1 

1 

CD  m 

1 

0) 

1 

1 

® X 

1 

•-  X 

C'* 

1 

1 

ro 

1 

1 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

•-  0 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

O 

o 

o 

o 

O 

CJ 

o 

O 

o 

o 

o 

o 

o 

o o 

o o 

o 

o 

o 

o 

O 

o o 

o o 

o 

0 0 

♦ 

♦ 

1 

X 

1 

X 

X 

X 

4 

X 

4 

X 

4 

1 

4 

4 

X 

X 

4 

4 

4 

X 

1 

X 

X 

X 

X 

X 

X 

X 

X 

1 

X 

X 

X X 

X- 

X 

X 

X X 

CD 

•- 

•- 

CT)  O 

9) 

rx  CO 

9 

rx 

•- 

ID 

9i 

o 

9) 

w 

X 

9 

rx 

o 

9 

9. 

o 

9. 

® 

9 

(N 

44 

CD  O 

o 

ro  X 

ro 

9 

cx 

0 0 

<7i  *I> 

«T» 

fji 

T 

ID  ro 

ro 

If)  If) 

ro 

9) 

9 

Cx 

■D  CO 

ro 

If) 

If) 

CD 

3. 

9 

rx 

CD 

X 

X if. 

in 

(V 

9. 

9 

4- 

X 

X in 

ro 

CD  9) 

ro  X o 

*0 

*- 

rx 

® 

ro 

If) 

•- 

o 

fx 

ro 

ro 

If. 

r- 

«- 

CJ 

•- 

(3> 

Cx 

ro 

ro 

if. 

1 

4- 

o 

4- 

0. 

Cx 

ro 

ro 

If) 

(- 

44 

o »- 

0) 

<X  X 

4- 

- a> 

r3 

cx 

If) 

o 

ro 

cx 

ro 

(O 

V- 

fO 

rx 

m 

o 

ro 

cx 

ro 

ro 

fO 

in 

o 

X 

rx 

X 

X 

4- 

X 

in 

o 

X 

rx  X 

ro 

X 

cx  CD 

cx 

r«. 

0* 

1 

1 

(XJ 

i/t 

r') 

1 

1 

1 

9) 

1 

t'- 

1 

tt) 

m 

ro 

1 

1 

1 

0) 

1 

1- 

1 

CD 

If) 

r> 

1 

1 

1 

9)  f" 

1 1 

ro 

in 

X 

1 

1 

1 

^ ro 

1 

o 

o 

o 

O 

o 

O 

O 

O 

O 

o 

O 

O 

o 

o 

o 

o 

O 

o 

o 

o — 

o 

o 

0 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

0 0 

♦ 

X' 

1 

4 

X 

X 

X 

4 

X 

X 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

X 

X 

X 

X 

X 

X 

X 

X 

4 

X 

X 

X 

X X 

X 

X 

X 

X X 

ro 

ir> 

rx 

m 

rj  (D 

X 

iT)  to 

O 

X 

•“ 

f». 

o 

r* 

ID  'O 

9> 

r*j 

X 

X 

4- 

9 

X If. 

ID 

9 

X 

ID 

cx 

X 

ro  9 

ID 

ro 

44 

T-  in 

ft) 

'D 

'0 

ro 

cr- 

If. 

■0 

ro 

3. 

r? 

t- 

rx 

CD 

if. 

ID 

ID 

r". 

o 

9 

1 • 

rx 

t ' 

if) 

If. 

'D 

ro 

rj 

9 

t' 

rx 

'7 

9 

ID 

X 

4- 

f'g  t- 

(\ 

9)  X 

CD 

ro  cx 

'0 

•- 

•o 

rx 

0. 

'1 

ro 

Cx 

'O 

ro 

0- 

f 4 

9 

3* 

9 

9) 

ro 

ro 

ro 

3. 

f X 

9 

3- 

9 

rn 

n 

ro 

X 

3. 

rx 

9 

9. 

9 

ro 

X 

X ro 

9. 

C'4 

9 

4- 

0 •- 

'I 

o 

•- 

r - 

ro 

o 

•- 

rj 

rx 

rx 

9 

if) 

f 

ro 

O 

CJ 

rx 

rx 

9 

If/ 

f >• 

X 

o 

r~ 

o 

rx 

rx 

9 

if. 

1 - 

X 

o 

4- 

O C'i 

(X 

9 

ro 

ro  cx 

rx 

r- 

(T\ 

ro 

X 

1 

1 

1 

f" 

! 

9) 

CO 

ID 

1 

1 

1 

X 

X 

CD 

,* 

i 

r-- 

1 

9> 

X (D  •- 

1 1 

1 

r* 

1 

in  •- 

1 1 

o 

o 

o 

O 

o 

o 

O 

O 

O 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

O 

o 

o 

o o 

o 

— 0 

o 

(J 

r> 

o 

r* 

r.) 

o 

o 

tj 

r_) 

rj 

rj 

o 

CJ 

o 

o 

CJ 

o 

CJ 

CJ 

r:j 

<J 

o 

rj 

o 

rj 

fJ 

fj 

CJ 

O 

o 

CJ 

o 

o o 

O 

o 

o 

0 0 

♦ 

X 

1 

4 

4 

4 

4 

4 

4 

* 

4 

4 

4 

X 

4 

4 

4 

X 

* 

4 

4 

4 

4 

* 

4 

4 

4 

4 

4 

4 

X 

4 

4 

X 4 

X 

X 

X 

X X 

rx 

O' 

r- 

» • 

'D 

ro 

n> 

r • 

X 

'D 

ro  9. 

t - 

LO 

O) 

If) 

rx 

9. 

f- 

X 

X 

X 

9. 

<x> 

'X) 

o 

in 

X 

9. 

ro  '0 

o 

9 9 

9) 

t'* 

ro 

h*  ro 

» • 

t 

f - 

ir. 

1- 

f > 

1 - 

0< 

'D 

1 - 

•r 

rx 

3» 

•D 

f - 

ro 

rx 

3. 

'D 

1- 

CC 

’T 

rx 

o 

‘D 

1 ' 

'D 

'7 

rx 

o 

‘D 

I'. 

ro 

9 

rx 

o 

•V  ( • 

ro 

9 

cx 

r- 

X 

CJ 

«- 

fj 

rx 

<:> 

•J 

1 • 

'll 

fO 

t. 

O 

f 

O. 

— 

ro 

'C; 

rj 

( 

00 

»- 

'0 

1' 

CJ 

J- 

O 

*- 

X 

1 

o 

l-^ 

i« 

ro 

1 • 

o 

ID 

4- 

ro 

4- 

9 rx 

rx 

«/) 

0. 

O 

»- 

rx 

ro 

X 

9' 

r,J 

•- 

CX 

ro 

9) 

O 

rx 

r^. 

X 

3. 

o 

4- 

rx 

X 

Tj 

3. 

o 

4> 

rx 

X ro 

9. 

o — 

CX 

X ro 

cx  0 

r) 

r- 

ro 

i 

r*) 

9* 

9)  ro 

X 

1 

9. 

1 

1 

1 

9) 

1 

X 

I 

ro  9) 

9) 

X 

3) 

1 

X 

1 

f 

1 

0) 

1 

X 

1 

X 

9) 

9) 

X X 9) 

1 I 

I 

1 

ro 

1 

— X 

o 

o 

O 

o 

O 

o 

O 

o 

O 

O 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

O 

4-  O 

o 

o o 

0 »- 

rj 

r> 

o 

rj 

O 

a 

rj 

o 

o 

O 

O 

c> 

rj 

O 

o 

rj 

o 

o 

r. 

r-j 

o 

O 

o 

CJ 

o 

o 

CJ 

o 

fj 

o 

o 

CJ 

o 

o o 

o 

o 

o 

0 0 

4 

» 

1 

1 

4 

4 

4 

4 

• 

4 

« 

♦ 

t 

• 

4 

4 

4 

4 

4 

4 

4 

X 

4 

4 

4 

4 

* 

4 

4 

4 

4 

* 

4 

X 

4 

4 ♦ 

4 

4 

X 

+ 4 

f-f 

'O 

t - 

'D 

r t 

'1, 

X 

ro 

", 

o 

9. 

'j 

9 

CO 

'D 

o 

f-' 

Cm 

CD 

If. 

x 

o 

cr 

if) 

ID 

rx 

ID 

Cj 

If. 

9 

rt 

*D 

3.  O 

o 

4- 

9 

ro  ro 

**y 

• 

•'j 

‘i. 

o 

f, 

o 

'i 

O’. 

fi 

", 

r"  . 

if. 

o 

o, 

'f, 

J. 

3i 

O 

■* 

if. 

'J 

if. 

3. 

r * , 

a. 

r. 

o 

if. 

3. 

ro 

9 Cj 

— 

if) 

ro 

'0 

- 9 

f J 

•- 

1 1 

' 4 

o- 

*0 

'‘4 

'f> 

X 

r ■» 

'J 

•1 

■/, 

• 

X 

'O 

O. 

0 

3. 

— 

0 

Ov 

4- 

'I' 

3. 

‘D 

4- 

■* 

f 

'1. 

ro 

»-  ' 

ro 

3. 

<D 

4- 

0 C J 

rx 

iD 

rx 

", 

•f. 

'Zi 

7. 

•- 

r , 

*0 

0. 

.. 

' X 

r*"/ 

0 

3- 

rx 

X 

CD 

4- 

ro 

'X 

*•. 

'*1 

4- 

r. 

X CD 

9. 

3. 

4- 

to  ro 

r) 

tv 

1 

• 

rj)  ro 

1 1 

CX 

X 

9)  ro 

cx 

1 

9) 

1 

1 

1 

9. 

1 

CO 

1 

rx 

X 

9) 

X 

CN 

1 

9> 

1 

,* 

1 

9i  X rx 

t 1 

9> 

»-  ro 

CO 

CX 

1 

ro 

1 

9 •* 

1 1 

o 

<•> 

o 

O 

o 

O 

o 

o 

o 

o 

O 

o 

o 

€) 

O 

o 

o 

o 

o 

O 

o 

o 

4-  0 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

O 

O 

o 

o 

o 

o 

O 

o 

o 

O 

O 

o 

O 

o 

o 

o 

O 

o o o 

o 

o 

o 

0 0 

♦ 

X 

1 

1 

4 

X 

4 

X 

4 

4 

4 

X 

4 

4 

4 

X 

4 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X- 

X 

X X 

X 

X 

X 

X X 

a« 

o 

t- 

‘D 

o. 

't 

ro 

in 

x 

r- 

rx  X 

ro 

ro 

r- 

f-- 

in 

X 

f'- 

•- 

X 

X 

X 

r- 

f * 

in 

9)  c- 

rx 

9»  9 

cx 

f'4 

r-  in 

ro 

c-*. 

cx 

X 9 

k/i 

o 

•f 

If. 

•T 

If. 

o 

i/. 

ro 

o 

(' 

'^1 

'D 

3i 

■V 

»- 

Cf 

r;, 

rj 

j ' 

3) 

'D 

3. 

ID 

•- 

X 

X 

o 

9. 

ID 

3. 

'D 

4- 

X 

ro 

o 

t « 

9. 

CD  c. 

'D 

4- 

ro 

ro 

ro  ID 

'jt 

T) 

rx 

•- 

ro 

r 4 

If. 

rx 

r*-) 

» - 

X 

•7 

9 

'I' 

( X 

If) 

rx 

ro 

r • 

'0 

'7 

9 

If; 

rx 

If. 

rx 

ro 

1 - 

a. 

9 

'7 

'D 

rx 

If) 

rx 

X 

r-  CD 

9 

9 

ID 

4- 

4-  0 

c* 

»<• 

l/» 

•- 

'0 

o 

ro 

O 

rx 

ro 

CX 

3J 

'0 

o 

ro 

o 

fX 

ro 

CX 

rn 

»- 

ID 

o 

X 

O 

CN 

X 

rx 

T) 

4- 

ID 

o 

X 

o cx 

X 

CX  CD 

4-  X 

t 

<o 

1 

1 

ro  r- 

X ID 

ro 

1 

r* 

1 

1 

1 

1 

1 

X ID 

1 1 

X 

r* 

9) 

ID 

X 

h* 

1 

1 1 

1 

1 

ro 

4-  X 

1 

o 

O 

O 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

o 

O 

o o 

O 

o 

o 

0 4- 

rj 

o 

ri 

o 

r_» 

rj 

o 

O 

CJ 

o 

O 

o 

o 

o 

o 

o 

o 

rj 

o 

o 

o 

o 

CJ 

o 

o 

o 

rj 

o 

CJ 

CJ 

<J 

O 

o 

o 

o 

CJ  CJ 

o 

o 

o 

0 0 

♦ 

X 

1 

1 

X 

4 

X 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

X 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

X 

X 

X 

X 

X 

X 4 

X 

X 

X 

X X 

O 

o 

rj 

h«. 

O 

ro 

o 

(7) 

•- 

X 

CD 

X 

rx 

o 

9. 

w 

rx 

o 

('• 

X 

o 

rx 

o 

0)  rx 

X 

f- 

X 

ro 

o 

rx 

o 

9) 

9 9 

9 

ro 

9 4- 

•- 

•7 

't 

-T 

lO 

»T 

O 

•7 

r • 

D 

rr) 

rT) 

o 

•- 

•- 

if) 

1 - 

D 

ro 

ro 

o 

9 

4- 

4- 

if) 

1 - 

ID 

ro 

ro 

o 

9 

O 

4- 

i/i 

I-' 

'D 

ro 

X 

O 9 

o 

4- 

ro 

9 

X 0 

CO 

rx 

•- 

ro 

'D 

ro 

♦T 

a* 

X 

’T 

rf 

ro 

rx 

'fJ 

•o 

-T 

9) 

(O 

’I 

9 

ro 

rx 

X 

ri 

9 

0* 

*7 

9 

ro 

rx 

CO 

ro 

9 

0) 

ro 

9 9 

fO 

CX 

ro 

4- 

CX  CX 

C- 

J) 

rx 

ro 

rx 

fj 

n 

If) 

o 

Cx 

ro 

rx 

O 

9 

4- 

m 

o 

I'. 

rx 

X 

fx 

o 

9 

4- 

.n 

o 

f'. 

fx 

X 

rx 

o 

9 •- 

ro 

o 

r* 

4-  X 

to 

'■ 

"■ 

c-* 

ro 

1 

ID 

1 

X 

1 

1 

1 

» 

1- 

1 

ro 

I 

ID 

X 

X 

1 

CD 

1 

9) 

1 

1 

1 

1 

r*  X 

1 1 

CD  ro 

n* 

1 1 

o 

o 

O 

o o o 

O 

o 

o 

O 

o 

o 

o 

o 

O 

o o 

o 

o o 

o 

o 

o 

o 

o 

o 

o o 

o 

O 

O 

O 

o 

o 

O 

o 

o 

o 

o 

o 

CJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o o 

o 

o 

o 

0 0 

X 

1 

1 

t 

X 

X 

4- 

4 

X 

4 

4 

4 

4 

4 

X 

X 

X 

X 

X 

X 

X 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

4 

X 

X 

X X 

X 

X 

X 

X 1 

<T> 

t' 

r- 

o 

ID 

rx 

X 

ro 

9) 

ID 

O ID 

ro 

fo  m 

(O 

m 

rx 

in 

rx 

T'- 

rx 

® 

4- 

X 

3) 

ID  O 

*o 

X X 

ro 

X in 

4- 

OJ 

ro 

4-  cx 

r>  iTi 

n 

ro 

ro 

rx 

ro 

•- 

ro 

o 

ro 

T. 

'D 

f'' 

9. 

O 

If. 

«- 

If) 

a> 

ro 

C) 

rx 

*• 

3. 

'D 

3- 

O' 

X 

9' 

■D 

r- 

3' 

O 

in 

4- 

If) 

9. 

X 

O' 

X — 

o>  '0  ro 

X ro  CO 

TO 

*- 

rx 

rx 

® 

ro  ffl 

•- 

ro 

If) 

'T 

a. 

r- 

rj. 

9 

9 

o 

ro 

C X 

If) 

ro 

9 

X 

rj 

ro 

if. 

4- 

9 

rr» 

4- 

9.  9 

9 

O 

X 

•- 

C-4  If. 

ro 

9 

X 

44 

9 9 

0» 

(O 

t- 

ID 

*• 

rx 

rf 

m 

ro 

o 

rx 

n 

X 

ro 

rx 

•- 

t'. 

9. 

3. 

t ' 

9 

4- 

rx 

9 

in 

X 

o 

rx 

rx 

X 

X 

X rx 

4- 

»-  9)  ro 

h- 

9 

4-  ID 

(T> 

o 

(T> 

ro 

m t' 

"■ 

X 

f'* 

m 

ro 

9)  X in 

1 1 1 

r* 

1 

0) 

1 

1 

1 

,■ 

1 

t 

i 

1 

1 

ro  n* 

1 1 

in  X 

1 1 

4-  X 

1 

O 

o 

O 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o — 

0 — 

O 

o 

o 

O 

O 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o o 

o 

o o 

o 

o 

O 

o 

o o 

o 

o 

o 

o 

o 

o 

o o o 

o o o 

0 0 

■f 

X 

1 

t 

X 

X 

X 

X 

X 

* 

4 

X 

X 

1 

4 

4 

X 

4 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

1 

X 

X 

X 

X X 

X 

X 

X 

X X 

o 

r* 

<J> 

lO  o 

ID 

X 

(T. 

o 

O 

X o 

X 

9 

•o 

-• 

•D 

•- 

rx 

m 

3i 

ID 

® 

o 

4- 

4- 

9. 

ro  9. 

X 

in 

00  9 

ID 

4- 

44 

h*  in 

rx  ‘O 

rx 

D 

fx 

rx 

rx  rx 

rx 

If. 

CT. 

rx 

X 

a> 

r- 

CO 

c. 

1- 

o 

9 

«- 

3. 

CJ 

•- 

‘D 

»• 

in 

a 

f4 

CJ 

r- 

X 

9. 

4- 

if. 

o 

9 

4-- 

0.  O 

4- 

<0 

«— 

cx  ID  CO 

Oi 

•- 

rx 

rx 

ro  CO 

n 

O 

ro 

•- 

o 

r^) 

rx 

Cx 

rx 

If) 

f' 

ro 

If. 

<7. 

V- 

3. 

9 

9 

o 

X 

»• 

4- 

X 

rx 

rx 

rx 

ro 

If. 

r-* 

X 

If)  0. 

44 

9. 

9 

4* 

ro  9 

r^ 

CD 

r» 

If) 

ro 

ro 

rx 

w 

c- 

X 

9)  f' 

9 

rx  cx 

9 

if)  X 

O 

rx 

rx 

X 

X 

X 

rx 

4- 

f'. 

3. 

0. 

I" 

9 

X 

rx 

9 

in 

X O 

(X 

rx 

X 

0 X 

a 

'■• 

'■• 

"• 

in 

ro 

9) 

t 

ro  in 

1 i 

t'* 

1 

X 

1 

1 

1 

1 

1 

1 

i 

1 

1 

O'  r~ 

1 1 

in  X 

1 i 

<7)  X m 

r**  ro  •- 

r.  4- 

1 1 

<y 

O 

o 

o 

o 

o 

o 

o o 

O 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o o o 

o o 

— 0 

o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o o 

o o 

o 

0 0 

« 

X 

» 

1 

X 

♦ 

♦ 

•4 

■* 

4 

4 

4 

X 

X 

4 

.4 

X 

4 

4 

X 

X 

4 

4 

4 

X 

X 

X 

♦ 

X 

4 

X 

X 

4 

X 

X 

X X 

4 

X 

X 

X X 

0) 

'J' 

a. 

r- 

*/» 

m 

m 

m 

iTi 

m 

lO 

IT 

lf> 

If. 

m 

in 

If) 

iT. 

iT. 

If. 

m 

III 

in 

in  in  tr) 

m 

If.  m 

m 

in 

in 

If. 

m 

in 

in 

in  in 

ro 

ro  ro 

X r* 

o 

■D 

rr 

II' d’ 

•D 

<D  H' 

o 

'D 

T" 

‘D 

"D 

T 

■D 

'D 

0 

’C 

•0 

■D 

D 

'X 

o 

'D 

‘D 

D 

'D 

'D 

CD 

'D 

ID 

ID  ID 

CD 

ID  (D 

9 

CD 

'0 

rx 

rx 

ro  go 

O 

ro 

fO 

ro 

ro 

n 

ro 

ro 

ro 

<0 

n 

ro 

ro 

ro 

lO 

ro 

n 

ro 

n 

X 

X 

X X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X X 

X 

X X 

44 

9 n* 

CJ* 

'0 

IT) 

in  If) 

lO 

If)  in 

fi 

X 

If. 

in 

If) 

If) 

.n 

If. 

.fl 

If) 

If) 

if) 

in 

•n 

in  in  in  in 

in  in 

1.0 

m 

in 

in 

■n 

in 

in  in 

ro 

ro  ro 

X — 

0> 

to 

o 

ui 

O) 

0>  (7)  9) 

XXX 

X 9) 

0) 

XXX 

X X X X X 

X X X X X X 

9>  9>  9)  9)  ® <T) 

ro  ro  ro 

ro  ro  ro 

0 

• M<n^ir»<or^ooa>o*-wfnntf>u)r^a)(y>o  — r'<n«T^n(ot^coa>o  — 

— •-•-C'»cNCKf'«C'<fM(Nr^c'<o*nr)r>tnr>r)cn 


r 


2^b2t1*00  -.6-J2i5*00  -.1I8?8‘0t  - . n!9^6---00  .13413-0'  .2'’C3-C1  .13486*01  .5-3031-00  -.12139  + 01  .32955-00 


(Ti  f- 
CJ  0)  if) 

n lo  f' 
04  O 

»-  O) 


- r>  fn  ( 
r)  O i - » 
O.  r - cj>  I 
o c£)  r\j  ' 
^ a»  »-  • 


• *7  'T  T 
I O 0i  «D 

• ">  o I- 
» CN  If)  'D 

> — t-  c« 


Cl  CD  a.  I 

O CD  fi*  ' 
^ If)  •- 


I CD  (N  (Ti 
• If)  l/» 


4 ♦ 4 ♦ 

Cl  ('■  n ^ 
C>  O O f') 
r)  CD  0. 
04  m CD  o 
t'  o»  - 


+ 44+44 


+ + 


+ + + + 


» 


♦ + 


I 


4 + 


I 


OCiTlDOVOrV'TOOf-*-  <7)Oa>r'-0>tD'Lf)f'i^O'7lf)J7)'7  ^ O t-  {D  Ot 
Lf<t'—ir.'-  — — t'-r4mt-'0(N4~i'(''CD'0»-'D('t-f^'D'D'Dr.i-'Do*r- 
CD  C)  CD  T X f*)  ID  CD  Oi  a<  O)  fo  Cl  CL  »J  'T  C i/»  (D  'D  T '7  (T.  ii  f-  CD  '7  1 0<  X X 

(Xj  04  rt  n CM  (\  in  ill  Ci  O)  t~>i  o cd  — c o •-  '^)  — i^.  7.  o — n — C'  o 


♦ ♦ ♦ + 

D>  (D  — ff) 

r-  ® ® ® 

X »j  *7  (J) 

— ® ^ r> 

^ ^ — ID 


I 


i I 


I 


I i 


o»-o-o- 
o o o o o < 


(^  O)  ^ O'  o* 
'7  ro  If)  r)  »r 


X If)  X 04  if) 
r)  '7  ® X X 

- X - 0)  - 
rs4  tn  04  CD  (N 

— If)  - U)  - 


+ ♦ 


4 + 4+  + + 44  + + 


+ 44 


XiX04XaiXXXh*'7X*7»-X®04Xf'4t'X«7XO  *-0 

X n X)  X '7  X X X <0  X X X t'  X X X X X X X X X X O)  X 04 
X — X •-  X — X »-  CD  •-  X — X •-  X — X »-  X X •-  X 0^  '7  X 

xo4rr)fNXxxo4Xxxxxxxxxo4XXfDr4m  xo 

J)  •-  U).  •-  If)  •-  If)  — ui  — X — X X — a)  +-  u)  *-  If)  «-  If)  »-  •- 


I 


00  — 
I o o o 

T + 4 

t - *7  X 
I-  X X 
— »7  X 
•“  O (N 
X X •- 


I 


I 


I I 


I 


+-  o o o 

O 0 J O O 


X «7  X *D 
04  '7  *7  If) 

a>  (■'4  if)  o 

04  T X X 

»-  r- 


*-  «-  o o 

O O O O ) 
4 + 44  + 


X *-  X X 
x X r * X 
*7  X X CO 
X o — o 
- ^ r)  «D 


4 4 4 4 

O O X X f" 

X X 7*  '7  X 

X 0<  X '7  O' 

04  0-4  'T  If)  f>+ 

- - f'  0» 


I 


•--00*- 

00000 

4 4 4 4 4 

X (D  X X X 

X X X X t-- 

•7  X X X X 

X O •-  0^  X 

— X X »- 


I I 


I 


I I 


I 


( i 


I I 


I 


I 


I t 


I 


I I 


I 


I 


I 


O Cj 

* 4 

X t • 
— 'T 
^ < < 
0^  7. 

--  X 


4 + ♦ + + + 


o»-o  — o--o*-o--o  — o»-o»-o*-o^  — o 

lOOOOOOOOOOOOOOOOOOOO 


+ 4 4 4 4 4 


4 + + + 


+ 444 


I 4 


oxcr)oxooi»-o»7aixx»-a)xa»ocj)xxxoxo)xa)xxr'-o  xt-'X'Dxf 

if.  ^ If)  --  CD  — ov  O'*  *-  •-  O •-  t'  •-  X »-  X *-  X *“  04  X 04  •-  O •-  X — O 04  X C'  •-  '4  t - X 1 

— <«  — r , ^ r4  »-  04  04  •-  04  — O •-  O •-  O'*  if)  00  — O f ' ( 

r-4  Ti  04  CT>  04  07'  X 01.  X 7i  X X X X X CTi  04  X X (J»  X X X (7i  X OTi  X oji  X oji  04  O O X X X •' 

— X»-X—  X--X*--X»-X^X»-X*-X*-X»-X»-X»-X»-X—  X»-  XX  — "“XC 


4 4 14  4 4 


1-  X X o> 
> 04  OJ  Gr  •- 
CD  t-  0 4 +- 
t 04  *-  0.2  X 
X X X +- 


X CD  CD 
X O O If. 
fj  t 'TO* 

04  X X *- 
— X X — 


4 4 1 4 4 + + 

»-  'O  X X X 04  O'- 

o X T <s^  — <0  in 

'T  I'  'D  f 4 — O 'O 
X X — ' O X X X 
•-  X X X r-  — X 


X ■ 


X *7  X 


"I  0>  ►-  0 Oj 

If.  7)  X t t ■ 

X 04  — O X 

^ ^ X — c. 


(D  X 0.  0. 
X »“  04  'D 
X 'I  04  ’7 
X 04  r2  X 
f-  — X -- 


X '7  X X 
X Gi  X 'O 
If)  D X t 
X 04  — O 
^ X - 


I I 


^7  X G) 

t - X *1 
X X 04 

o.  r-  »- 


i I 


I t 


I 


I I 


t I 


I 


I 


I 


I I 


I 


I I 


t 


X G)  X (0)  X <0  f - X cc  X G.  C7i  If.  X X 'O  X '7  X oo 

'7  'D  04  G.  O 'CJ  O.  CD  - 'll  'D  — G.  O X 0>  X t-  'O  'X)  iT) 

04  -r  J)  CD  X t - x X G 04  ; jj  x c-  x x 04'7X+- 

O0GX,04—  OXXXOX'DX—  O'nx  oo-t- 

x+-»-*-x»-o.  r'+-x+~^»“X  — Oir^  •-Gi»- 


I O 02  « 
4 » 

X f' 

'7  X 
OG  D 

04 


G.  -D  *7 
G.  'T  If. 
G.  G)  — 
— — X 
If)  X ^ 


I-  'I  ~ 

O X Gi 
CD  Oi 


X G - '■4 
t>  X G.  ‘7 
fG  D 02  *7 
X 04  O 04 
»-  h-  G - 


X X O X 
'D  t- 


4 4 4 + 

O )D  (T)  X 
X X G X 
If)  C2  X X . 
— f+.  ^ 04 
G.  ^ ^ 


I I 


t I 


1 I 


I I 


I 


I I 


I 


».»*♦  + ♦♦!  »»44»»+4|  * r 

Sj  ■>  if)  O f - ^ ’ 0^  D G 'i.'  Oj  f..  - X X G O if.  — 

— ",  G.  :•  'X>  f - »->  '7  G — X X ''  'G  f ' X ’T  o^  r, 

• f • » f < if,  O O'*  •/.  G - ' «L>  04  l/l  0.2  X 1'4  rn  '* 

X G - f.  X - 'D  "j  C2  X G - If.  c:  — X X O X ^0 

^G-®®--G*-0'l--G'-®ro—  G^x»-G 


+ + 4 + 4I  •t44-+»4 

• 0 4 X C)  0 . X •-  t - X G>  07.  f)  o 

•7  0 4 0:  - X G I '•  o V.  ' 

- fv  if.  O 0)  ‘ > X • . - 0 * u-  I 

if)X  — X.XOX/G-X  . — 

•XX+-G  — ON  — G — ax  — 


I I 


Gi  G G.  G 
: "<  r - G.. 


I 


1 I 


lT.  X XX  0 4 

X 0>  — O'-  X — 


‘♦♦♦♦■♦•♦♦♦I  ♦♦4* 

. c:  o r-»  f - +-  " . x X X Gi  i - 04  ( - X 

X " G X X 0 O'.  J".  X G J<  X t- 

G u.  '•  X 'D  ' . ■ - : ' . ' . 0.  X 

•X-X—  '-.0*.  XOG'‘X04Gi 

04X  — G.  — — — xx  — — xx  — r-x 


I I 


0’  O.  0 4 o>  r - 

o — oot- 

f - o X X Gi 
G,  G X Xi  04 
If)  C»  - O)  — 


I I 


G.  « • O C 
OD  iD  X r 

r-  — 1 


I I 


I I 


+ + + + + + 


+ + 


I 


I 


» 


I ( 


I 


I 


{ I 


+ + 


+ + 


+ + + + 


c.  X I 

X 0-4  • 
I 'D  lO  { 

G)  X 1 


) X X O'-  — (D  — h»  Gi  If)  O X G)  f"  O O)  G.  f"-  G (D  04  O X 04 

) if.  I-  — O — O G.  I-  Gi  O f-  X G G)  X If.  — X)  ’ifM  - G 04  X X 0-4  X 

. — O G t-  O X X G.  X X O X — '0  04  G ('  G>  04  OD  if.  Gj  X X X O') 

— X ")  X G 04  X 04  X X X OO  0*4  (D  X X X O 04  X O 04  X O 04 

*-x— XX  — — — — xx  — X— 


. O O O O O ) 

+ + + + 

— '7  X G»  X 
X 04  X X 04 
If.  O)  X X X 
X O OV  X)  C) 


+ + + + + 

— — X — oo 

X X 04  X X 

oo  If)  G)  X X 

04  X O 04  X) 


I 


+ + + + + 

ON  0»  oD  X O oj) 

04  t-  X ON  X 

X X G.  X X 

cf  c\  in  o 


I t 


X X If)  G» 
X — ( • 
CD  oj  O')  G 
— O O OG 
^ ^ X - 


I I t 


i I I 


I 


I 


I 


C)  04  rj 
O - X 
Gi  G)  *7 
— 0:2  -• 


.000 

4 4 4 

G - I- 

Gi  X X 
t - fG  if) 
f*.  X X 
04  - X 


o - - • 
Oj  o o « 

+ + I 
0.2  If.  04 
Ifi  — X 
O X If) 
O 04  G 
G)  — G 


+ + 
r>  ('  f-  r-~  I 
ON  04  G G) 
X 04  04  — 
04  OG  — X 

— m r->  — 


• o — o — — o—  00  — — — o 
)OOO0j  000  02  00000 

♦ + + 4 


+ 4 


4 + + + 


XOXXXXONX  — GXXaC'—  GIOONXG  — XiGiX  — GXOXxX04(DXf'O0Na) 
04  G.  *7  X I ' '7  O X X X.  X — — G —XGX0NO00  — G.'0  — r'X)04r'-t-C4X'Dr4X0G  — OX  — 

O X O 04  G o.  'D  G t--  X X O If)  o X 04  X X G X O')  '7  t - X G t'  G '7  t - '7  G I ' '7  G t ' G G 0-+  G G 

G.  — OOX  — O — ('XXOONGOnX  XXXiDXXiDXXXXXXXXinXXXXXXX 

— — G — — ON  — XG)  — G — x X — r-'X  — t'-x  — t-'-x  — r-if)  — r'-x  — r''if)»--o-'*x  — 


I I I 


I I 


I 


I 


I 


I 


I 


I 


I 


I 


000000000000000000000000000000000 


G X)  X 

^ X ^ 

— X — 


X X G' 
G-  G If) 
G If.  7) 


X O X 
Gi  G.  X 


^ X - 
^ X — 


X X Gi  X 
G.  G.  X X 
G If.  X G 

— X - — 

— X - — 


I 


+ ♦ • 


I 


I 


+ J 4 ♦ J + + I 


+ + +I  44+4-4  + 


+ + ♦ + + + 


— GXXG*  X — GXX®X  — Xr4--Xr4GX  GXt-^XX  — GXCNGXOOXOGXOGXXOX 
X X 0.  G.  If)  a.  t - ifi  G'  O'  If)  O'  X X G.  O*  X Gi  CN  If)  X X — CN  X — t-'  G X r4  X X •-  X — X a t-  G fN  — C4  t-  — 

X X G If)  X 71  X X G X X G X X '7  If)  X '7  X X f4  X X X C.  '7  (-*  G If.  Gi  X X X X X X G X — X X r4  X G 

X — — X—  — X—  — X—  — X—  —X—  — «D—  — OX(TiGrNG)r4XGXX—  OXXXO—  OXXG 

X—  — X—  — X—  — X—  —X  — — X — — X—  »-f4  — XX  — CN  — XX  — CN—  — G—  G — — CN  — X® 


I I 


I I 


i I 


I I I 


0-00  — 00-00  — 00  — 00-00  — 00-00  — 00-00- 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOl 


♦ + 


+ + + + + 4 + + 


4 ♦ ♦ 4 4 + + 


r#—  XOr4XGXX7'  %*•  — — GXr*^.  XXGXXXt-»  — — XCNT-if) 

O G>  — G.  Gy  X X X X ('  X X r-  X X X X X.  X X if)  X X X If.  X X G X (D  X 

X G r*  X G X G f4  If.  G rs  if.  G r4  if.  G c , x G X x G X if.  G f 4 X G r4  If. 

OXGOXGOXGOXGOXGOXGOXG  OXG  OXGOX'7  O 
r - xt-^  — Xf'  — — XT--  — ®r-  — XT'*  — Xh-  — Xf*- 


— G X 
X t*  f4  - 

G X rj  x 

X G G' 
X X 


I 


o — o — — o 
000000 

+ 4 + + 4 + 


I f 


I I 


I 


I I 


f4  X 'G 
X X O X o 
X t'  X — X 
X O — O X 
— X — — G 


— o — — o 

O O O O O ( 

+ + + + + + + + I 


— cr  Gi  X X 

X — 0»  G — 

G X G If.  O 

X X O — X 

— G — — CN 


X G X X O 
r4  — o r-  X 
ro  — f-  X X 
X X X fN  X 
— X X — CN 


— 00  — 0 — — 
0000000 
+ + + +♦  + ♦ 
OXGX  — X — 
X X X X fN  — X 
».  Xt-Xh*®- 
fV  G ® X O — O 
X X — fN  — — 


I I I 


--00  — 
o o o o o 

4 4 4 4 + 

G X X O T' 
G.  G r-*  G r- 
a.  - X X 
— 4*4  O CN  X 

^ ^ ff  ^ 


o — — 
000 

4 I 4 
G G +• 
X G.  r4 
— r*  t - 
4*.  X O 

o.  X — 


I I 


— 00  — — 
00000 

♦ 4 4 4 + 

t - f - X X o 

J-.  r.  f - X 

— ® O ® X 

X — X X — 

— X X - — 


I I I 


I I I 


00  — 

000 


o o — — o 
o o o o o 


X O TN 
X G r- 

- G G 
I-  — X 

- ® - 


♦ + 444  + + + 


0-0-- 
O O O O ( 


^xxx  — 

G 'N  X -7  X 
f'+XC-G  — 

— a>  — o 

® fN  — — G 


I 


X f-*  m X 
X t'-  X — 
X X — fv 
X X G>  X 
h-  •-  (J)  X 


X X 
I"*  — 
O X 


I I 


00- 

>000 

4 + 4 

X x 

m - r- 

X O X 
— fN  X 
X If)  - 


I 


I 1 


I 


I I 


I 


I 


-00 

000 

4 4 + 

^ t-  — 

X X G 
X — G 

- - ® 


o o 

4 I 
X <’2 

X G G 
X G 
0) 


— o — 00  — o — — — o — 00  — o — — — o — 
ooooooooooooooooooooo 


I 


4 4 


4 4 


I + + 


XOG®®Gr4XrNXfN®Gi{T»G  — G»®GG® 
xor^fNXxx  — xxot-  — (Nxx  — xa  — r-'* 

OG-t'f4—  — Glf)OOi  — r'+X  — — GXO®  — 
XG  — XG)r4f''XGXG  — X(T.X^•XGXG  — 
— lO—  OT**—  X—  ®—  X—  X— 


I 


I I 


I I 


I t 


1 I 


I 


I I 


XG  X®f«®®0’ 


fNXGX®f^*®®0—  CNXGX®f-®®0  — CNXGX® 
— — — — ——  — — CNCNrNfNCNfNfNfNfNCNXXXXXXr) 


•-fNXGX®r-^®®o- 


rNXGX®f«-®®0  — (MX 

— — — — — — — — CMXCMn 


B-22 


J 


557 


o t » o o o o o c>  o o o o 

l*»4**l  «**♦♦« 

O'  ^ r*>  I ■ '£«  f - rt>  o»  O)  1/1  fO  o 

-tot  o»  ’ 4 * 0>  r J O)  m f - 'D  Oi 

; '1;  a<  'T  '7  o-  r cd  *7  •»  — 

O.  » * »•  n •-  r*)  f-i  - ro  •-  <n  •- 

fO  ^ ^ 't  f-  ^ <0  ^ 


III  I I I I 


««******»44f4^ 

CTi  »-  ftj  'C*  f'*  t fO  r,j  fO  t/j  r«  JO 

'T  *7  <i>  'D  ’*  *7  f ^ r - fO  'O  •-  in 

O.  "I  *7  ’U  *7  f4  fO  0>  *7  7' 

r'4  ■?  If.  f - r*’*  f>  •■•  rj  CM  f4  *7  i/)  t>- 

. (7.  - f'>  (O  » - 0> 

I I I I I I I 


***«l  * * * * * ^ ¥ ^ 

(tj  rj  «D  'U  i/>  f i/»  • ’1  'Xi  f'i  04 

fiUi'  . o.  *7  i/>  7;  'O  lO  Oi  in 

•f  r I : - 'i  J.  ’1  t 4 ('i  'D 

•r  ^ f 4 ' > *■ ) <o  f 4 "t  'ft  r) 

-oi^oto—  - mfo  — --a* 


I 


I 


I I 


I 


•-  'f)  -■  1-  it)  *7  f*  o.  a<  rj  a,  fo  «o 

••  'o  f 4 ' • 'n  < , ti,  •j  itj  -t  fo 

- n ► I -'  *.»:•>  -7  oi  I 'D 

f 4 • 

l/j  '0  - - 'O  o.  - r • *7  ^ C4  (t> 


I 


I 


I 


I 


I 


' 4 ' 4 - "I  CJ 

^ iti  to  — 


ir»  If)  -0  o.  iri  'u  t ' 0'  o.  *•  > rx.  - - tr.  ••  *7  ••o  r 4 rj  u,  ) . rj,  at  oj  m oi  •/  'Xi  •- 

't  ' . ift  '/  f j ■',  — I 'i;  f . f')  ■:  if.  If  iX.  ■£;  O.  f'*  If.  »■  > t -7  - O,  »■  > fO  ‘I-  ' - 'X) 

f 4 ^ O'  trjrjfjt^fjfjtj4,^^c^r4'‘,-r’*  ; f''  X-  0'  o 

",  ; . ' ; . f.  »/J  •.'<  i/»  in  '•••,/.  ,f.  . : i/.  . . .r.  ,>./<./>  if.  ./  i/»  if-  i/»  if»  iT.  ifi  .f » ''4 

O'  o.  O’!  O'  ot  (7i  0>  (7)  0»  O'  O'  O'  O'  O'  O'  O'  O)  O'  O'  O'  O'  O*  O*  0»  O'  O'  O'  O'  01  O'  O'  Oi  O'  O'  o>  Oi 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


I 


• o c»  o < 

♦ ♦ f 

I O I/)  - r 

'X.  'Ti  '4  ' 

> '*’>  (X.  'h  • 

rj  n i[,  fj  r4  .x.  o ‘‘4  ifj  o < 


O '■>  - O '*>  O O ^ ^ ^ ^ ^ ^ O O O O r-  o o o o 

oocj  0000  00000000  OOOOOOOOOf? 


, ^ o, 

, r 4 'f.  'Tj  f I ■ 


I o)  in  *7 

•J  it  't* 


♦ 1 ♦ 4 


4ff44  + 44  ♦♦♦I 


4 1 ♦ 4 f 4 4 ♦ 


0000000 


'!>  O'  'D  — 03  — no  fo  •-  o n f - n o O'  in  — cd  n o o O'  n f ••  n to  'o 

fo  m '7  f ’t  '0  in  n *7  o.  ift  o '7  n.  »d  f * o-  ’7  c:,>  n O'  — t - o ( 4 ix.  'O  < 4 

n '7  — n O'  o o iX'  'ft  • ' 'O  ixi  tx,  ^ 'jj  ’i  a t iX'  o n 1 4 *7  » - t - 

fO  o o 'o  'fj  o o.  o — r4  n n n x^  *-  ex  o ;fi  »7  — - ix’i  rn  0i  n o •-  r4 

»T04fnx'<'7  7«aD«-'-»-^*-^  — "---tDiDooxeNXN’rtnco---  — 


444 
01  xQ  in 

O'  <0X4  7- 

X4  -7  r4  in 
n xn  n 04 


Ci 

1 

1 

0 

0 

1 

0 

0 

1 

CJ 

0 

0 

1 

0 

1 

1 

i 

7 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

r » 

c-> 

X> 

Oj 

o» 

r t 

X J 

r ) 

( f 

XJ 

X > 

0 

0 

C ) 

XJ 

X> 

0 

XJ 

0 

0 

0 

Xj 

CJ 

OJ 

XJ 

0 

0 

CJ 

0> 

CJ 

xj 

0 

0 

0 

0 

0 

OJ 

OJ 

XJ 

CJ 

OJ 

0 

0 

0 J 

CJ 

0 

( 

* 

4 

4 

t 

4 

4 

4 

4 

4 

4 

4 

4 

« 

4 

4 

4 

4 

4 

4 

f 

4 

4 

4 

4 

4 

1 

1 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

» 

4 

4 

O' 

m 

O' 

fX> 

XT) 

o> 

n 

'T 

O' 

xn 

r • 

XT) 

fl; 

XXJ 

If 

n 

•0 

n 

n 

r 4 

IX» 

lX) 

0 

ON 

0 

ON 

ID 

fii 

XJ 

04 

n 

0 

■7 

04 

•- 

XJ 

(D 

'7 

CJ 

7- 

CD 

6 

8 

5 

O' 

•0  n 

-7 

'X, 

•7 

xn 

•7 

'0 

X4 

iX' 

'7 

X4 

iXi 

'I . 

J' 

O' 

IXI 

iXi 

*7 

'7 

'C> 

'JJ 

( 4 

7 - 

'O 

•- 

04 

'0 

n 

0. 

o> 

XT' 

x:. 

ft) 

7 

f ■ 

XJ 

X4 

•7 

1' 

7 • 

•- 

I/' 

ifi 

•f, 

iX' 

'7 

‘D 

rX' 

'7 

'i> 

iX' 

*7 

0) 

tXi 

'7 

'0 

n 

Oi 

•- 

■D 

O' 

x 

XX) 

'7 

r- 

'7 

fX) 

0 

'f 

to 

0 

•- 

♦7 

'O 

O' 

'JJ 

'7 

r ' 

O' 

xfj 

n 

'7 

-7 

'O 

'7 

CJ 

•- 

X*) 

n 

'7 

•7 

T' 

n 

in 

n 

'X. 

iX> 

n 

Oi 

iX. 

n 

O) 

r4 

_ 

C> 

X 4 

Xf 

n 

n 

f4 

0. 

0 

1 - 

n 

XT' 

r- 

04 

•- 

rj 

ON 

n 

n 

n 

X 4 

04 

X4 

0 

1' 

•• 

O' 

r~ 

iXi 

1 

r 

m 

1 

r- 

i/i 

1 

f ♦ 

m 

t- 

1 - 

1 

O' 

i 

1 

1 

1 

O' 

1 

<D 

i 

*7 

1 

'O 

1 

04 

m 

1 

e-' 

O' 

1 

1 

1 

1 

1 

O) 

CD 

1 

70 

•7 

1 

CO 

ON 

1 

in 

0 r> 

0 

0 

rj 

0 

0 

XJ 

0 

0 

OJ 

0 

0 

0 

0 

X) 

7r> 

0 

0 

0 

r> 

n 

r> 

0 

rj 

0 

0 

0 

0 

x^ 

r> 

0 

0 

0 

O) 

0 

0 

0 

0 

0 OJ 

0 

0 

0 

0 

0 

0 

0 

rj 

0 

CJ 

r? 

fj 

0 

0 

0 

0 

0 

0 

OJ 

OJ 

CJ 

0 

0 

CJ 

0 

0 

0 

» 

1 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

♦ 

4 

4 

4 

« 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4- 

04  XD 

Oj 

n 

0 

0 

'T 

0 

O' 

r 

n 

04 

•- 

iXl 

04 

7- 

*- 

w 

Oi 

n 

m 

V- 

«o 

70 

0 (D 

•- 

0 

03 

in 

*- 

in 

7' 

70 

XT) 

70 

OJ 

(D 

XD 

m 

XD 

0 

CD 

»7 

ON 

'0 

r 

n 

r - 

n 

'n 

» 

n 

•7 

Oi 

03 

'7 

n 10 

CJ 

l/i 

n 

Xj 

i/»  'D 

7- 

XU  'li 

•- 

0 

1 ' 

04 

If 

0 

1 

xjj 

1' 

04 

03 

X,J 

*7 

OJ 

0 

71) 

O' 

ru 

'7 

7 - 

O' 

O' 

i/'  0 

•7 

70 

r . 

O' 

/' 

f > 

n 

tX. 

If 

'7 

O' 

n 

r - 

iXi 

r 

c> 

'O 

«- 

If) 

'7 

1 • 

If 

O' 

04 

m 

ff) 

If) 

-7 

I - 

IX* 

O' 

ixj 

1X( 

O' 

•- 

IXI 

’7 

r • 

iXi 

r4 

'D  CJ 

•7 

•- 

(X) 

*7 

04 

0. 

r4 

'O 

'7 

") 

•> 

0 

n 

X'l 

XJD  0 

04 

0 

o> 

'0 

0 

m 

0 

04 

n 

X'4 

0 

O'  <0 

0 

'7 

0 

04 

n 

04 

0 O' 

'O 

Oj 

0 

0 

04 

n 

04 

0 

0 

(O 

0 

'7 

OJ 

ON 

n 

*■ 

04 

1 

1 

a> 

•0 

1 

C4 

1 

0 

1 

»7 

1 

1 

( 

*0 

1 

04 

04 

(X) 

1 

1 

1 

<n 

1 

r4 

04 

(O 

1 

1 

1 

(O 

1 

ON 

ON 

(O 

1 

t 

1 

O' 

\ 

ON 

ON 

CO 

1 

1 

0 

0 

0 

0 

0 

0 

x:^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

XJ 

0 

0 

0 

0 

rj 

0 

OJ 

0 

0 

X3 

0 

0 

0 

0 

0 

r> 

0 

0 

D 

0 

0 

0 

0 

0 

0 

0 0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 0 

0 

0 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

1 

•a  n 

rf 

•D 

X; 

O) 

O' 

•- 

«. 

n 

n 

•7 

IX 1 

0 

n 

Xf 

<0  0- 

O' 

w~ 

O' 

r» 

0 

ID  xf 

-7 

0 

0 

10 

0 

m 

1/1 

XO 

•7 

773 

•7 

m 

Oi 

m 

m 

on 

ON 

O) 

t- 

•- 

Xj 

Tj 

n 

0 

*7 

0 

n 

'7 

04 

T> 

’f 

f- 

0. 

04 

-7 

0. 

0,# 

O' 

o-j 

O' 

0 

w- 

0. 

•- 

'li 

•> 

0 

•- 

OJ 

03 

X4 

r- 

(/' 

1 - 

0 

0. 

0 

O' 

CX* 

m 

04 

n 

0 

•“ 

0 

*- 

to 

•- 

r 

'7 

iX' 

i/. 

iX' 

»- 

x^' 

' J 

O'  iX. 

0- 

1 

If 

xx) 

'O 

i/* 

X4 

Xj 

'7 

0 

n 

'1 

'7 

'L» 

10 

04 

'7 

•- 

O'; 

'7 

•- 

10 

ID  04 

CJ 

♦7 

0 

X4 

n 

X4 

7 - 

7' 

74 

0 

•- 

n 

7" 

n 

•n 

n 

in 

Oi 

0 4 

f • 

04 

0 

r . 

0. 

n 

n 

74 

n 

0 

O' 

7- 

0. 

«- 

n 

fO 

04 

n 

-7 

O' 

t- 

0* 

•- 

n 

n 

74 

Xf 

•7 

XD 

7- 

(D 

*- 

xn 

n 

ON 

ro  *7 

•7 

1 

0 

1 

1 

04 

1 

to  (S 

1 

04 

t 

1 

n 

r* 

1 

7 

1 

O' 

in  m 

1 

n 

1 

1 

1 

O)  in  ® 

1 

n 

1 

1 

1 

Oi  in 

CO 

n t' 

1 

1 

1 

CD 

m 

1 

O) 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

r> 

0 

x-> 

0 

r) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Oj 

0 

0 

0 

0 

0 

0 

0 0 

r> 

03 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

000 

0 

♦ 

4 

4 

1 

4 

4 

4 

4 

1 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

4 

4 

1 

4 

4 

r> 

n 

'D 

rr 

0 

•n 

XT. 

O) 

0 

Oi 

X4 

04 

0 

04 

0 

— 

»-• 

•D  <0  *0 

ri 

T» 

x:> 

n i--- 

Xf 

'D  TO 

ON 

m 

xf 

•- 

Xf 

0 

•- 

74 

04 

Xf 

CD 

O' 

CD 

CD 

r-i 

X* 

O' 

It 

f 

•7 

X4 

0 

rr 

0 

0 

'0 

'X> 

T 

O' 

7 • 

to  04 

lf»  O' 

'/i 

0 

» - 

O' 

If 

7 - 

7. 

f'. 

7' 

Uj 

*7 

7- 

0* 

70 

'X> 

■7 

1' 

’■1 

•“ 

(O  ® n 

'O 

r 

04 

. 

0 

'XI 

0 

O' 

•• 

r 

X4 

XO 

n 

1 

X4 

n 

O' 

n 

•J 

O' 

n 

•7. 

0 

n 

*7 

XN 

O' 

0 

Xf 

iXi 

n 

•- 

O' 

0 

n 

If) 

n 

•- 

O' 

•7 

n 

to 

7- 

0 

O' 

•7 

rr, 

7 • 

xn 

n 

•4 

XT' 

X4 

r- 

X4 

r- 

XO 

04 

Ifl  04 

n 

04 

'O 

04 

n 

04 

If 

XN 

xf 

04 

If) 

X4 

n 

04 

*7 

XN 

(Ti 

n 

tn 

fO 

® m K 

70  in 

"x; 

(O  m 

r* 

CO  in 

h* 

CO 

.5 

.7 

CD  m 

4 4 4 4 4 


•-  4 or  4 


•yint04'®®»o*-fN4^»TU><o 

f^f^cNXNO^oinDfnooxnror) 


c - 


B-23 


COI  E*T-NAME  UNT  EC  QOr  SEC  C0,.CC  NEXT  lIVIT  READ  *RITTEN 


1 ACTIVE  DEVICES  ( 0 FULL) 

TP-OPS,  6 MRITES.  1 READS.  10488  WORDS  XFD 


APPENDIX  C 

USER  INFORMATION  FOR  THE  AUGMENTED  MATRIX  PREPROCESSOR  AUCMAT 


The  following  Includes  modifications  which  were  made  under  separate  contract  to 
the  Naval  Surface  Weapons  Center,  Contract  Number  N60921-77-C-01 1 2 . 

This  appendix  includes  a copy  of  the  users  manual,  and  a sample  input  deck  and 
subsequent  output  for  the  Infinite  cylindrical  shell  problem  presented  In  Section  4. 


C-1 


THIS  functional  component  of  the  UNCER/JATER  shock  analysis  code 
CCEPTS  data  from  ThE  FLUID  MASS  MATRIX  PROCESSOR  AND  THE 
S’’RUCTURAL  ANALYZER  TO  CONSTRUCT  THE  SPECIFIC  CONSTANTS  AND  ARRAYS 
That  are  used  in  the  staggered  solution  procedure  for  the 
TRANSIENT  Response  analysis  of  submerged  structures 
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rtic’  followint'  <1 1 sciifiK  Ion  If)  provldfd  ns  nn  aid  to  nsor  understanding  of  the  sample 
output  that  Is  inelnded  here. 

After  a summary  of  the  fluid  mesh  geometry  arrays  (see  Appenriix  li)  the  first  Item 
needing  explanation  is  that  entitled  "Fluid  Klement  Wetted  Freedom  Indleator".  This  Is 
siniiily  a listing  of  the  Input  variable  NDIFOS  (see  user  manual)  for  ea<h  fluid  element. 

ITie  section  "Struetural  Crid  Point  Numbers  Associated  With  Internal  Sefpience  Num- 
bers" contains  a correspondence  table  that  relates  the  Internal  serpjence  numbers  assigned 
by  the  fluid  mass  processor  with  the  external  structural  node  number  assigned  by  the 
user . 

The  next  Item  entitled  "firld  Point  and  Freedom  Number  for  Fach  Row  of  .Stiffness 
Matrix"  identifies  an  Integer  vector  tliat  Is  constructed  by  the  user  In  the  Skyline 
Utility  (see  Figure  1-1,  also  Appendix  F).  For  each  structural  equation  the  entry  In 
the  vector  consists  of  ten  times  the  structural  node  number  plus  the  local  degree  of 
freedom  number. 

The  last  Item  requiring  explanation  is  the  "Freedom/Kquat Ion  Correspondence  Table". 
Tills  Is  an  Integer  matrix  of  6 rows  and  as  many  columns  as  there  arc  structural  node 
points.  Any  particular  row  corresponds  to  a local  degree  of  freedom  number  while  a col- 
umn corres(>onds  to  the  Internal  sequence  number  for  a particular  external  node  number. 

The  matrix  entry  for  any  particular  set  of  row  and  column  is  the  structural  equation 
ntimher  for  that  pair. 
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SUMMARY  OF  DATA  STORED  ON  PERMANENT  FILE 
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APPKNDIX  D 


USKR  INFORMATION  FOR  THE  TIME  INTEGRATION  PROCESSOR  TIMINT 


The  following  includes  modifications  which  were  made  under  separate  contract  to 
the  Naval  Surface  Weapons  Center,  Contract  Number  N60921-77-C-01 1 2. 

This  appendix  Includes  a copy  of  the  users  manual,  and  a sample  input  deck  and 
subsequent  output  for  the  infinite  cylindrical  shell  problem  presented  In  Section  4. 
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UNIT  OPERATION  CONSTANT  IN  SECONDS,  CONSISTING  OF 


FLOATING  addition.  A FLOATING  VULTIPLY.  AND  INDEYIN' 
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POST  PROCESSING  (SUBROUTINE  RESDSP): 


The  following  discussion  is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  included  here. 

First,  the  amount  of  storage  required  for  the  run  given  in  the  output  refers  solely 
to  the  blank  common  that  is  set  in  the  main  program,  UNWASH.  An  error  exit  is  taken  if 
insufficient  storage  is  available  and  the  user  must  see  that  more  is  provided  either  by 
a recompilation  on  UNIVAC  llOO-OS  or  by  a field  length  request  on  CDC. 

Sector  address  information  for  the  response  and  restart  files  that  is  listed  at 
various  places  in  the  output  is  extremely  important  for  subsequent  restart  runs. 

Although  transient  response  results  can  be  displayed  as  part  of  the  run  such  output 
has  been  deferred  to  the  post  processing  phase  in  Appendix  E for  this  sample  problem. 
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RESTART  DATA  FOR  T • 3.000000  WRITTEN  AT  LOCATION  121  ON  PERMANENT  FILE  CYL'RESTAR 
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APPENDIX  E 

USER  INFORMATION  FOR  THE  POSTPROCESSOR  POSTPR 

The  following  includes  modifications  which  were  made  under  separate  contract  to 
the  Naval  Surface  Weapons  Center,  Contract  Number  N60921-77-C-0112 . 

This  appendix  includes  a copy  of  the  users  manual,  and  a sample  input  deck  and 
subsequent  output  for  the  inflnte  cylindrical  shell  problem  presented  in  Section  4. 
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this  ru\CTioNAL  component  of  the  u'.derwater  Shock  analysis  code 

IS  RESPONSIBLE  FOR  THE  TABULATION  AND  PRInTER-PlOT  GRAPHIC  DISPLAY 
OF  SELECTED  TRANSIENT  RESPONSES  AND  PSEUDO- VE LOC 1 T Y SHOCK  SPECTRA 
UPON  COMPLETION  OF  AN  UNDtR*lATER  SHOCK  ANALYSIS  OF  A SUBMERGED 
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NWETHS  I NUMBER  OF  STRUCTURAL  HISTORIES  (EITHER 

DISPLACEMENTS  OR  VELOCITIES)  TO  BE 
DISPLAYED  FOR  WHICH  THE  APPROPRIATE 
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FACTOR  E.F  VALUE  OF  MULTIPLICATIVE  LENGTH  CONVERSION 

factor  to  be  applieo  ■^0  The  displayed 
structural  transient  response  histories 


SHLIST  L TRoE  if  PSEuDO-VELOClTy  ShCCK  SPECTRA  ARE 

TO  nE  listed  in  tabular  form,  otherwise 
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The  following  discussion  Is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  included  here. 

Perhaps  the  only  item  needing  discussion  is  the  transient  response  tabular 
listings.  The  desired  responses  are  displayed  In  matrix  form  so  that  each  row  contains 

the  entire  history  of  a particular  degree  of  freedom  except  for  the  first  row  which  is  ' 

time.  Each  column  therefore  contains  the  instantaneous  values  of  the  complete  set  of 
response  variables  desired  at  a particular  time.  Each  row  Is  identified  by  the  struc- 
tural or  fluid  node  and  its  degree  of  freedom.  The  letters  D,  V,  and  P stand  for  dis- 
placement, velocity,  and  pressure,  respectively.  > 

Essentially  the  same  format  is  used  for  the  pseudo-velocity  shock  spectra  except 
that  the  first  row  is  now  frequency  rather  than  time. 
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TRANSIENT  RESPONSE  HISTORIES: 


1/0  P . 48338*00  .48843+00  .49380+00  .49876+00  .50277+00  .50515+00  .50686+00  .50900+00  .51160+00  .51401+00 
0/0  p .68044+00  .63747+00  .73945+00  .65715+00  .61509+00  .59371+00  .58605+00  .59083+00  .60688+00  .63137+00 
9/0  P .20016+00  .27911+00  .34677+00  .39782+00  .43517+00  .46655+00  .49693+00  .52523+00  .54825+00  .56570+00 
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APPENDIX  F 


USER  INSTRUCTIONS  FOR  INTERFACING  WITH  USA 

To  use  tS.e  Underwater  Shock  Analysis  (USA)  Code  in  its  linear  stand-alone  mode,  the 
user  must  first  construct  a permanent  data  file  that  contains  the  structural  mass  and 
stiffness  matrices  and  some  assorted  bookkeeping  information.  The  purpose  of  this  appen- 
dix is  to  describe  the  structure  of  the  file  and  to  specify  how  it  is  to  be  created.  At 
this  time  utility  routines  that  carry  out  this  task  have  been  written  for  SPAR  and  NAS- 
TRAN.  An  abbreviated  form  of  this  file  is  also  required  when  USA  is  coupled  with  a non- 
linear structural  analyzer  and  such  an  interface  also  exists  for  STAGS. 

USA  contains  the  data  management  utility  module  DMGASP  that  carries  out  all  data 
transfer  activities  between  core  and  peripheral  storage.  This  is  done  by  unformatted 
and  unbuffered  data  transmissions  and  it  is  imperative  that  DMGASP  be  used  to  create  the 
structural  interface  file.  Otherwise  the  user  must  supply  or  have  access  to  a similar 
means  of  direct  transfer.  Section  3 of  [10]  contains  a comprehensive  discussion  of  the 
half-dozen  or  so  DMGASP  commands  that  are  required  to  activate,  position,  write  upon, 
read  from,  and  free  a peripheral  storage  device.  Subsidiary  commands  also  exist  for 
error  handling  and  listing  of  selected  information  pertaining  to  auxiliary  storage. 

The  current  configuration  of  USA  uses  a diagonal  mass  matrix  associated  with  a 
lumped  mass  representation  of  the  structure,  assumes  that  there  is  no  velocity  dependent 
structural  damping,  and  further,  only,  single  precision  matrices  may  be  processed.  In 
addition,  if  the  stiffness  matrix  has  been  reordered  or  reduced  in  any  way  for  input 
to  USA  the  mass  matrix  must  also  be  reordered  or  condensed  so  that  its  degrees  of 
freedom  (DOF)  are  the  same  and  appear  in  the  same  order  as  in  the  stiffenss  matrix. 
Finally,  the  stiffness  matrix  must  be  placed  in  a multi-block*  skyline  format  as  dis- 
cussed in  [13].  This  description  consists  of  a Matrix  Master  Record  (MMR)  followed  by 
a series  of  Matrix  Value  Records  (MVR)  which  contain  the  numerical  values  of  the  matrix. 
These  are  the  only  constructs  the  user  need  be  concerned  with;  all  others  required  are 
already  embedded  in  USA.  During  construction  of  the  MMR  a logical  device  index  (LDI) 
must  be  set  in  the  record  which  USA  will  access  later.  For  UNIVAC  operation,  this  should 
be  set  equal  to  twenty  (20),  while  for  GDC  operation  this  should  be  set  as  two  (2). 

The  file  structure  required  is  shown  in  Table  F-1  where  NDOF  stands  for  the  number 
of  structural  DOF  which  USA  must  process.  NMMR  is  the  number  of  words  in  the  matrix 
master  record,  and  NWBL  is  the  number  of  words  in  each  matrix  value  record  (which  is  ex- 
pected to  be  the  same  for  each  record).  NWBL  should  also  be  an  Integer  multiple  of  448 
for  most  efficient  use  of  storage. 

* For  small  problems  a single  block  is  permissible. 


F-1 


Number  of  Words 

Data 

1 

NDOF 

NDOF 

Diagonal  Mass  Matrix 

1 

NDOF 

NDOF 

Grid  Point/DOF  Vector 

NMMR 

NMMR 

Matrix  Master  Record  for 
Stiffness  Matrix 

NWBL 

First  Matrix  Value  Record 
for  Stiffness  Matrix 

NWBL 

Second  Matrix  Value  Record 

for  Stiffness  Matrix 

NWBL 

Last  Matrix  Value  Record 
for  Stiffness  Matrix 

The  Grid  Polnt/DOF  vector  consists  of  an  Integer  value  for  each  global  DOF  from  1 
through  NDOF  that  is  constructed  as  ten  times  the  external  node  number  plus  the  local 
DOF  number  that  apply  to  that  particular  structural  equation. 

For  example,  if  the  87th  DOF  to  appear  in  the  mass  and  stiffness  matrices  corres- 
ponds to  the  second  degree  of  freedom  at  a node  identified  externally  as  4637  then  the 
87th  entry  in  the  Grid  Polnt/DOF  vector  would  be  46372.  Local  translational  degrees  of 
freedom  should  be  numbered  1-3,  rotational  degrees  of  freedom  should  be  numbered  from 
4-6  and  any  others  should  be  numbered  with  7-9.  If  more  than  9 degrees  of  freedom  are 
carried  at  any  node  it  is  a simple  matter  to  change  the  factor  of  ten  to  one  hundred  in 
a few  places  in  USA  to  accommodate  this. 

It  should  be  noted  that  records  1-4  are  accessed  by  the  USA  pre-processor  AUGMM 
before  the  time  integration  phase  of  the  analysis  commences.  This  portion  of  the  file 
is  required  for  both  USA  in  the  linear  stand-alone  mode,  and  for  USA  when  it  is  inter- 
faced with  a nonlinear  structural  analyzer.  In  this  latter  case,  the  fifth  and  succeed- 
ing records  do  not  exist. 
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